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Abstract. We consider the 3d cubic focusing nonlinear Schrodinger equation 
(NLS) idtu + Au+ \u\'^u = 0, which appears as a model in condensed matter theory 
and plasma physics. We construct a family of axially symmetric solutions, corre- 
sponding to an open set in -ff^xiaiC^^) initial data, that blow-up in finite time with 
singular set a circle in xy plane. Our construction is modeled on Raphael's construc- 
tion [33^ of a family of solutions to the 2d quintic focusing NLS, idtU+Au+\u\'^u = 0, 
that blow-up on a circle. 



1. Introduction 
Consider, in dimension n > 1, the p-power focusing nonlinear Schrodinger 

(1.1) idtu + Au + = . 
This equation obeys the scahng symmetry 

u{x,t) solves (11. ip =^ \~u{\x,\^t) solves (11 .ip . 

which implies that the homogeneous Sobolev norm is scale invariant provided 
s = f — ^jzj- The equation (II. ip also obeys mass, energy, and momentum conservation, 
which are respectively defined as 

11 r 

M[u] = \\u\\l2 , E[u] = -\\Vu\\l2 - ^^\\u\\l^l, , P[u] = Im J Vuudx. 

In the subcritical setting (1 < p < 1 -|- ■:^), there exist soliton solutions u{x, t) = 
e**Q(x), where 

(1.2) -Q + AQ+\Qr'Q = 0. 

We take Q to be the unique, radial, positive smooth solution (in M") of this nonlinear 
elliptic equation of minimal mass. For further properties see, for example, |14j-[T5] 
and references therein. 

The local theory in {p < 1 + ;^) is known from work of Ginibre-Velo [10]. Local 
existence in time extends to the maximal interval {—T^, T*), and if T* or are finite, 
it is said that the corresponding solution blows up in finite time. The existence of blow 
up solutions are known and the history goes back to work of Vlasov-Petrishev-Talanov 
'71 [IQ], Zakharov '72 |12] and Glassey '77 |l2j who showed that negative energy solu- 
tions with finite variance, ||xtio||L2 < oo, blow up in finite time. Ogawa-Tsutsumi [31] 
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extended this result to radial solutions by localizing the variance. Martel [22] showed 
that further relaxation to the nonisotropic finite variance or radiality only in some of 
the variables guarantees that negative energy solutions blow up in finite time. First 
result for nonradial infinite variance (negative energy) solutions blowing up in finite 
time was by Glangetas - Merle [llj using a concentration-compactness method (see 
also Nawa [30] for a similar result). Positive energy blow up solutions are also known 
and go back to [IQ], [12] (see [381 Theorem 5.1] for the precise statement). Turitsyn 
[39] and Kuznetsov et al [20] extended the blow up criteria for finite variance solu- 
tions of subcritical NLS and, in particular, for the 3d cubic NLS showed that finite 
variance solutions blow up in finite time provided they are under the 'mass-energy' 
threshold M[mo]-E'[mo] < c(||Q||i2), where Q is the ground state solution of (11. 2p . and 
further assumption on the initial size of the gradient ||Vmo||l2 > c{M[uq], WQWl^)- 
Independently, the same conditions for finite variance as well as for radial data (for 
all subcritical, supercritical NLS) were obtained in [13]. [Similar sufficient 
blow up conditions for finite variance or radial data for the critical NLS (s = 1 or 
p = 1 + are due to Kenig - Merle [18]; for the situation in the critical NLS 
(s = or p = 1 + ^) with data refer to Weinstein [H] for the sharp threshold and 
Merle [22] for the characterization of minimal mass blow up solutions.] The nonradial 
infinite variance solutions of the 3d cubic NLS blow up in finite or infinite time pro- 
vided they are under the 'mass-energy' threshold having the same condition on the 
size of the gradient as discussed above. This was shown using a variance of the rigidity 
/concentration-compactness method in [16], thus, extending the result of Glangetas- 
Merle to positive energy solutions. Further extensions of sufficient conditions for 
finite time blow up were done by Lushnikov [21] and Holmer-Platte-Roudenko [13] 
which include blow up solutions above the 'mass-energy' threshold and are given via 
variance (or its localized version) and its first derivative if the data is not real. It is 
also possible to construct solutions which blow up in one time direction and scatter 
or approach the soliton solution (up to the symmetries of the equation) in the other 
time direction, see [13], [5], [6]. 

A detailed description of the dynamics of blow-up solutions in the critical [p = 
1 + ^) case has been developed by Merle-Raphael [25l [28l [231 [271 [261 [32] . They show 
that blow-up solutions are, to leading order in H^, described by the profile Qb{t), 
rescaled at rate A(t) ~ ((T - t)/log | log(T - ^)|)^/^ where b{t) ~ 1/log | log(T - t)\. 
Here, Qt is a slight modification of Q - see notational item N[5] in ^ for the details 
of the definition in the 2d cubic case. 

In the supercritical, subcritical regime (l-|-^<p<l + ;^), one has a 
large family of blow-up solutions as discussed above, but there are fewer results char- 
acterizing the dynamical behavior of blow-up solutions - we are only aware of three: 
Raphael [33j (quintic NLS in 2d) and the extension by Raphael-Szeftel [33] (quin- 
tic NLS in all dimensions) and Merle- Raphael-Szeftel [29] (slighly mass-supercritical 
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NLS). In this papeio, we consider the 3d cubic equation {n = 3, p = 3) which is 
a physically relevant case in condensed matter and plasma physics. We adapt the 
method introduced by Raphael [33] to give a construction of a family of finite-time 
blow-up solutions that blow-up on a circle in the xy plane. Raphael constructed a 
family of finite-time blow-up solutions to the 2d quintic equation {n = 2, p = 5) which 
is if^/^ critical, that blow-up on a circle. He accomplished this by introducing a radial 
symmetry assumption, and at the focal point of blow-up (without loss r = 1), the 
equation effectively reduces to the Id quintic NLS, which is critical and for which 
there is a well- developed theory characterizing the dynamics of blow-up. We employ 
a similar dimensional reduction scheme - starting from the 3d cubic NLS (which is 
H^^^ critical) we impose an axial symmetry assumption, and construct blow-up so- 
lutions with a focal point at {r,z) = (1,0), where the equation in (r, z) coordinates 
effectively reduces to the 2d cubic NLS equation, which is critical. 
Our main result is 

Theorem 1.1. Let Qb = Qb{r,z) be as defined in i\0 in ^ There exists an open 
set V in if^xiail^'^)' defined precisely in ^ such that the following holds true. Let 
Uq E V. Then the corresponding solution u{t) to 3d cubic NLS ( (11.11) with n = 3, 
p = 3) blows-up in finite time < T < oo according to the following dynamics. 

(1) Description of the singularity formation. There exists X{t) > 0, r{t) > 0, z{t), 
and 7(t) G M such that, if we define 




1 



Qb 



r — r{t) z — z{t) 



m ' m 



) 



u{t) 



Uit) - Mcore(t) , 



then 



(1.3) 



u(t) — u 



m L^{W^) ast^T, 



(1.4) 



VM(t)||L2(K3) 



Vncorc(i)||L2(]R3) ■ 



1 



log(T-t)|'= 



and the position of the singular circle converges: 



(1.5) 
(1.6) 



r{t) r{T) >0 ast^T, 
z{t) z{T) ast^T. 



A result similar to the one presented in this paper was simultaneously developed by Zwiers , 
see remarks at the end of the introduction. 
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(2) Estimate on the blow-up speed. We have, as t T , 

f T-t \ 
'''' ^Wn .og|.og(T-.)| ) • 

(1.9) 7(«)~|log(T-()|log|log(T-«)|. 

(3) Structure of the remainder u* . For all R> small enough, 

1 



(1.10) / \u*ir)\'drdzr^- 

J|r-r{T)|2+|^-2(T)|2<i?2 l^Og I ^OgK\j 

and, in particular, u* ^ forp > 2. 
(4) iJ^/^ gain of regularity outside the singular circle. For any R> 0, 

(1.11) u* e H^/Wr - r(T)|2 + |z - ;2(r)|2 > . 

A key ingredient in exploiting the cylindrical geometry away from the z-axis is 
the axially symmetric Gagliardo-Nirenberg inequality, which we prove in ^ This 
takes the role of the radial Gagliardo-Nirenberg inequality of Strauss [36] employed 
by Raphael. In ^ we collect most of the notation employed, and in ^ we outline 
the structure of the proof of Theorem 11.11 Most of the argument is a lengthy boot- 
strap, and in ^ we enumerate the bootstrap input statements (BSI 1-8) and the 
corresponding bootstrap output statements (BSO 1-8). As the output statements 
are stronger than the input statements, we conclude that all the BSO assertions hold 
for the full time interval of existence. The steps involved in deducing BSO 1-8 under 
the assumptions BSI 1-8 are outlined in ^ and carried out in detail in ^^4211 ^|2T] 
already proves (11. lip in Theorem II. H the if gain of regularity outside the singular 
circle. The proof of Theorem II. II is completed in ^ }22] - [Ml In ^22\ we prove the log-log 
rate of blow-up (II. 7p . In §231 we prove the convergence of the position of the singular 
circle, (II. 5p and (II. 6p . The proof of the size estimates on the remainder profile, (ll.lOp 
is the same as for Theorem 3 in [26] (we refer the reader to Sections 3, 4, 5 there). 
Finally, in §241 we prove the convergence of the remainder in L^, the estimate (II. 3p . 

A similar result to Theorem 11.11 but for a slightly smaller class of initial data was 
recently obtained by Zwiers [43] (our i^axiai class of initial data is replaced by a 
smoother version H^^:^^^). We mention that the methods to prove our main theorem 
do not treat the cubic equation in higher dimensions, but it is addressed by Zwiers in 
[i3| Theorem 1.3]. We think that it should be possible to adapt our method to treat 
the critical (p=l + 4/(n-l)) case in higher dimensions, proving the existence 
of solutions blowing up on a circle (dimension 1 blow-up set), however, dealing with 
fractional nonlinearities is a delicate matter. Zwiers' result treats blow up sets of 
codimension 2 in all dimensions; our approach might be able to treat the blow up 
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sets of dimension 1 in all dimensions. (In 3 dimensions, of course, codimension 2 
equals dimension 1, and thus, the results intersect.) 

1.1. Acknowledgements. J.H. is partially supported by a Sloan fellowship and NSF 
grant DMS-0901582. S.R. is partially supported by NSF grant DMS-0808081. This 
project has started at the MSRI program "Nonlinear Dispersive Equations" in Fall 
2005 and both authors are grateful for the support and stimulating workshops during 
that program. We thank Ian Zwiers and Fabrice Planchon for remarks on this paper. 

2. A Gagliardo-Nirenberg inequality for axially symmetric functions 

We begin with an axially symmetric Gagliardo-Nirenberg estimate, analogous to 
the radial Gagliardo-Nirenberg estimate of Strauss [36]. Consider a function / = 
f{x,y,z) = f{rcos9,rsm9,z) = f{r,z), independent of 9. We call such a function 
axially symmetric. 

Lemma 2.1. Suppose that f is axially symmetric. Then for each e > 0, 

(2-1) ll/lli».)<;ll/llii,.(.>.)liv/iii.^^(.>.). 

Proof. The proof is modeled on the classical proof of the Sobolev estimates. We 
use the notation V = {dx,dy,dz) (i.e., not {dr,dz)). Note that \\f{x,y,z)\\L2^^ = 
||/(r, ^)||l2^ ^ , and also note drf{r, z) = {dxf){r, z) cos 6' + {dyf){r, z) sin6', and thus, 
\drf\ < |V/|. Observe that for a fixed r > e, z G M, by the fundamental theorem of 
calculus and the Cauchy-Schwarz inequality, 

|/(r,;.)r = |/(+oo,z)p-|/(r,z)|2 

r+oo 

<2 / \f{r\z)\\drf{r\z)\dr 

Jr'=e 

oo 



<- / \f{r\z)\\drf{r\z)\r'dr' 



e 



<-||/(^',^)llL^,,(.'>e)l|V/(r',;.)|L. 

^ r' dr' ^ ' r' dr' ^ ' 

and also, 

/• + 00 

l/(r,^)P<2 / |/(r,/)PJ(r,;^')M^'< ll/(r,^OIlL^J|V/(r,z')llL^,- 

J z'=—oo 

By multiplying the above two inequalities, and then integrating against rdrdz, we get 
[ \f{r,z)\'rdrdz <-( [ ||/(r', V/(r', rfz' 

Jr>e,z ^ \J z ^ / 

(r+oo 
J \\fir,z')\\,.jVfir,z')\\,ydr 

Following through with Cauchy-Schwarz in each of the two integrals gives (12.11) . □ 
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As a corollary of (12. ip . we have for 1 < p < 3, that 



op-l 

(2-2) ll/llSt.„>.,=^Tirll/llk.<.->.>ll^/llS.(.>.)- 

This follows by the interpolation estimate (Holder's inequality) 

ll/llSi<ll/ll£l/llS"'- 

Before proceeding, we present a simple application of Lemma I2.1[ 

Corollary 2.2. If u is a cylindrically symmetric solution to 

idtu + Am + \u\^~^u = 

forp < 3 m that blows-up at finite time T > 0, then blow-up must occur along the 
z-axis. Specifically for any fixed e > 0, 

(2.3) hm||VM|U2^^(,<,) = +00. 

Proof. Fix any e > 0. All norms will be with respect to dxdydz. For any t > 0, 
by ( 12. 2p . we have for u = u{t) 

P ~r ^ e 2 

Using the inequahty a/3 < ^^a^^ + we obtain 

II v^lli. < + -^lklli^i.(.<.) + 4rll«lll^ • 

Since Ymit/'+oo \\^u{t)\\L2 = +oo, we obtain that lim^/i+oo ||'"(^) ||Lp+i(?->e) — By 
the (standard) Gagliardo-Nirenberg inequality, we obtain (12. 3p . □ 

3. Notation 

Recall that we will impose the axial symmetry assumption, i.e., 

u{r, 6, z) = u{r cos 6, r sin 6, z) 

is assumed independent of 9. The equation (II. ip in cylindrical coordinates (r, z), 
assuming the axial symmetry, is 

(3.1) idtU + -dru + dlu + dlu + \u\^ m = 0. 

r 

Denote by Q = Q(?"), t = x G M^, a ground state solution to the 2d nonlinear 
elliptic equation (which corresponds to the mass-critical cubic NLS equation in 2d): 

-Q + AK2Q + |Qpg = 0. 

We emphasize that Q is a two-dimensional object. 
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We now enumerate our notational conventions. 

N 1. We will adopt the convention that V = {dx,dy,dz) is the full gradient and 
A = d^ + dy + d"^ is the full Laplacian. When viewing an axially symmetry function as 
a function of r and z, we will write the corresponding operators as V(^r,z) = {dr, dz) and 
A.[r,z) = dl + dl- Note that under the axial symmetry assumption, A = /S.(r,z)+'r^^dr. 

N 2. In the argument, parameters A(t), 7(t), r{t) and z{t) emerge. Rescaled time, 
for given h^, is 



(3.2) s{t) = I T^ + So, So = e«o 



A2(t') 



and rescaled position is 



-/ N f — r(t) _ , , z — z(t) 

We note that the r and z can both be negative, although there is the restriction that 
X{t)r + r{t) > 0. Introduce the full (rescaled) radial variable 



We have r = X{t)r + r{t). To help avoid confusion between r and f, we will use the 
following notation: given two parameters A(t) and r{t), define 

/iA(t),r(t)(r) = (A(t)f + r(^))l{A(^)f+r(t)>o}• 
Often the X(t),r{t) subscript will be dropped. 

N 3. The inner product (■, ■) will mean 2d real inner product in the f, z variables. 

N 4. If / = /(f, z), then define A/ = / + (f, z) ■ [df, dz)f, the generator for scaling. 
Observe that Af = f + Rd^f. 

N 5. For a parameter b, Qb{r,z) is the 2d localized self- similar profile. Specifically, 
following Prop. 8 in [2^ which is a refinement of Prop. 1 in |28], given b, rj > 0, 
define 

R, = ^^VT^- 

There exist universal constants C > 0, rj* > such that the following holds true. For 
alio < rj < rj*, there exists constants e*(?7) > 0, b*{ri) > going to zero as ?7 — )• such 
that for all |6| < b*{ri), there exists a unique radm/ solution Qb (i.e., Qb depends only 
on R) to 

f A^f,z)Qb -Qb + ibAQb + QbW = 

QbiR)e^ > for G [0,Rb) 

Qb{Rb) = 
[Q(0)-e*{r])<Qb{0)<Q{0) + e*{r]). 
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N 6. For a parameter b, Qb{r,z) = Qb{R) is a truncation of the 2d locahzed self- 
similar profile Qh- Specifically, (following Prop. 8 in [27j), given 6 > 0, > small, 
let = y^l — rjRf, so that R^ < R^. Let (pb be a radial smooth cut-off function such 
that (pbix) = for \x\ > Rb and 0fe(x) = 1 for |x| < R^ , and everywhere < 0(x) < 1, 
such that 

Now set 



()'Jloo + II A^bllioo ^ as |6| ^ . 
Qb{R) = Qb{R)UR) ■ 



Then 

(3.3) A(f,5)Qf> -Qb + ibAQb + Qb\Qb\^ = -^6, 
where 

(3.4) - ^6 = A(pb Qb + 2VQb ■ V(pb + ibQb (r, z) ■ V(pb + (vl - Vb)Qb\Qb 
with the property that for any polynomial P{y) and A; = 0, 1 

(3.5) ||P(z/)vI/f^lU^ <e-^-/l^l. 
In terms of \&;,(f, 5), we define an adjusted \^;,(t,f, 5) as 

(3.6) ^b{t, r, z) = ^b{t, r, z) ^^^^dfQb{f, z) 



so that Qb solves 



(3.7) A(^f,z)Qb + -dfQb -Qb + ibAQb + QblOb]"^ = -^b. 

/i 

We also split Qb into real and imaginary parts as 

Qb = s + ie. 

It is implicitly understood that S and G depend on b (or b{t)), and when we want 
to emphasize dependence, it will be stated explicitly; this decomposition is done only 
for the truncated profile Qb (not for Qb). Similarly, we denote \E'fe = Re\E' + iIm\E' 
and = Re + i Im ^f. 

N 7. The Qb satisfies the following properties. 

QP 1. ((44) in [28j) Uniform closeness to the ground state. For a fixed universal 
constant C > 0, 

||e^^(4-Q)llc3 ->0 as6->0. 

In particular, this implies that ||e'^'^(9^B||i2(-f:^^) — as 6 — ?■ for < A; < 3. 

QP 2. ((45) in [28j) Uniform closeness of the derivative dbQb to the ground state. 
For a fixed universal constant C > 0, 



e^^'idbQb + tlR'Q) 



as 6 -)> . 

C2 
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QP 3. (Prop. 2(ii) in [28j and Prop, l(iii) in [2l]) Degeneracy of the energy and 
momentum. Specifically, 

(3.8) 1^(0^)1 < e-(i+^'')^/l^ since 2E{Q,) = j A^^Ofc, 
and 

(3.9) Im J V(^f,z)Qb 4 = 0, Im y (f, ~z) ■ Vif,z)Qb 4 = -^\\RQb\\h. 

QP 4. The profile Qf, has supercritical mass, and more precisely 

d 



< 



d{b^ 



fe2=0 



j IQfeP = do< +00. 



QP 5. Algebraic relations corresponding to Galilean, conformal and scaling invari- 
ances: 

(3.10) A{RQb) - RQb + tbRAQb + RQb\Qb\^ = 2df,Qb - Mb 

(3.11) A(i?2Q,) - i?2g^ + ^hFe^Q^, + feQb\Qb? = 4:AQb - 

(3.12) A(Agfe) - AQb + iAQb)\Qb\^ + 2Qb(S(AS) + 0(Ae)) 

= 2{Qb - ibAQb - ^b) - A^b - ibA^Qb. 

The proof of the above identities are similar to Prop. 2(iii) in [28] adapted to the 
2d case. For example, to obtain the third equation from scaling invariance, multiply 
(13. 7p by and take argument to be (Af, Xz), then differentiate with respect to A and 
evaluate the derivative at A = 1. Note that ^|a=i (A^^6(Af, Xz)) = 3^b+(f, z)-V^b = 
2\E'fe + A^'ft, and thus, we obtain the claimed equation. 

N 8. The linear operator close to Qb in (13.71) is M = (M+, M_), where (e = ei + ^£2) 

A / 2S^ \ 

M+(e) = -A(,,,)ei - -9,ei + d - + 1 |4rei - 2See2, 

\\Qb\^ J 

X ( 2B^ \ 

M_(e) = -A(,,,)e2 - -d.t^ + £2 - --^ + 1 \Qb?t2 - 2me,. 

WQbr / 

N 9. For a given parameter b, the function (b{'^, z) is the 2d linear outgoing radiation. 
Specifically, following Lemma 15 in [27], there exist universal constants C > and 
?7* > such that for all < < r^*, there exists 6* (77) > such that for all < 6 < 
b*{vi)^ the following holds true: There exists a unique radial solution Cb to 

{d} + dl)Cb-Cb + ibACb = ^b 



j \{df,di)Cb\'^ dfdz < +00, 
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where is the error in the Qb equation above. The number Tb is defined as the 
radiative asymptotic, i.e., 

Tb= Jim \R\''\Cb{R)\'. 

The Cf)(^^; and Tb have the following properties: 
ZP 1. Control (and hence, smallness by ZP2) of (b in 

|2 ^ -pl-Cr; 



y iV(f,5)C6 
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ZP 2. Smallness of the radiative asymptotic 

y\R\ > Rb, e-(^+^'')"/' < -r, < \R\\biR)\^ < e-(l-^'')"/^ 

5 



ESETIEIEEI. We note 



N 10. We will make the spectral assumption made in [25 
that it involves Q and not Qb- A numerically assisted prooflis given in Fibich-Merle- 
Raphael |9]. Let (see N 4) 

AQ = Q + (f, z) ■ VQ, A'Q = AQ + (r, 5) ■ V(AQ). 

Recall that (-, ■) denotes the 2d inner product in (f, z). Consider the two Schrodinger 
operators 

Ci = -A + 3Q[(f,z) ■ VQ], 
C2 = ~A + Q[{r,z)-VQ], 
and the real- valued quadratic form for e = ei + ie2 E H^: 

H{e, e) = (£iei, ei) + (£262, £2) . 

Then there exists a universal constant 61 > such that Ve G -ff^, if the following 
orthogonality conditions hold: 

(ei,Q) = 0, (ei,AQ) = 0, (ei,yQ) = 0, 
(£2, AQ) = 0, (£2, A^Q) = 0, (e2, VQ) = 0, 

then 

H{e,e)>~6i (^j \Ve\^ + j |e|Vl^l^ . 
N 11. The ring cutoffs are the following. The tight external cutoff is 



1 for < r < yI and for r > 
for i| < r < i|. 

16 — — 15 



The spectral property was proved in Id and in dimensions 2-4 has a numericany assisted proof. 
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1 for < r < T and for r > 4 
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^'^""^ ^0 fori<r<2. 



The internal cutoff' is 



,0 for < r < 7 and for r > 3 
r for I < r < 2. 



4. Description of the set V of initial data 

We now write down the assumptions on the initial data set V. Let V be the set of 
axiaUy symmetric Uq G H^{E?) of the form 

^*o(r, z) = —Qbo — , — ;^ — e ^° + UQ[r, z) , 

AO \ AO AO / 

and define the rescaled error 

eo(f, z) ^ Ao'Uo(Ao(f, 5) + (ro, 2;o))e"*^° , 

with the following controls: 

IDA 1. Localization of the singular circle: 

|ro — 1| < a* , l^ol < ■ 

IDA 2. Smallness of bg, or closeness of Qh^ to Q on the singular circle: 

< 6o < a* . 

IDA 3. Orthogonality conditions on eo: 

(Reeo, \R\^^bo) + (Imeo, |^|'06„) = , 
(Reeo, {r,z)T,bo) + (Imeo, {r,z)Qbo) = 0, 
-(Re eo, A'OfeJ + (Im eo, A'EftJ = , 
-(Re eo, AObJ + (Im eo, AEfeJ = . 
IDA 4. Smallness condition on cq: 

£(0) = [ |V(,,,)eo|VAo,ro(r) cifdi + / |eo(f, ijl^e-'^l cifd^ < lJ, . 

-'l(?,z)l<10/feo 

IDA 5. Normalization of the energy and localized momentum (recall ip from N 11): 



AqI^'oI + Ao 



Im / Vip • Vuq Uq 
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IDA 6. Log-log regime: 

0<Ao<exp(-exp(|^ 
IDA 7. Global smallness: 

||mo||l2 < a* • 
IDA 8. i/^/^ smallness outside the singular circle: 

\\XoUo\\h^/2 < a* . 

IDA 9. Closeness to the 2d ground state mass: 
IDA 10. Negative energy: 

E{uo) < 0. 

IDA 11. Axial symmetry of uq. 
Lemma 4.1. The set V is nonempty. 

Proof. This follows as in Remark 3 of Raphael [33] . □ 

5. Outline 

Now let u{t) be the solution to NLS with initial data from the above set P, and let 
T > be its maximal time of existence (which at this point could be +cxd). Because 
Mo is axially symmetric and the Laplacian is rotationally invariant in the xy plane, 
the solution u{t) will be axially symmetric. Thus, we occasionally write u{r, t). The 
first step is to obtain a "geometrical description" of the solution. Since we are not 
truly in the critical setting and cannot appeal to the variational characterization of 
Q, we need to incorporate this geometrical description into the bootstrap argument. 
The lemma that we need, which follows from the implicit function theorem, is: 

Lemma 5.1 (cf. Merle-Raphael |2B] Lemma 2, [2S] Lemma 1). If for < t < ti, 

there exist parameters {X{t),^{t),f{t),z(t),b(t)) such that ||e(t)||//i < a^, ^ 1 on 
<t <ti, where 

e(f, z) = e*^A'u(Af + f,\z + z) — Qi{r, z), 
then there exist modified parameters {\(t),'y(t),r(t), z(t),b(t)) such that e defined by 

e(f, z) = e*'^Au(Af + r, Xz + z) — Qb{r, z) 
satisfies the following orthogonality conditions (with e = ei + ie2): 
ORTH 1. (ei, + (e2, W^m) = 

ORTH 2. (ei, (f , 5)S6(t)) + (es, (r, ~z)Qb(t)) = 
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ORTH 3. -(ei, A2e,(i)) + (e2,A2S,(i)) = 
ORTH 4. -(ei, AQb{t)) + (e2, A^bit)) = and 

\Vu{t)\\L2 



livgiu^ 

where 6{ao) — )■ as ao 0. 



+ \b{t)\<6{ao), 



Note that the condition ORTH 2 is a vector equation. 
We also use the notation 

It is important that we consider, by default, u and u as 3d objects in the spatial vari- 
ables. Thus, when we write \\^L2,'ffe mean \\uJ^x,y,z,t)\\L2(dxdydz) = \\u{r, z,t)\\L2{rdrdz)- 
On the other hand, we consider, by default, Qfe(t) and e to be 2d objects in the spatial 
variables, and thus, if we were to write ||Qb(t)||L2, we would just mean \\Qb{t){f^, z)\\L2(df^dz)■ 
}^oweveT, if working with the r, z variables and the function Qb{t) and e, we will write 
the integrals out to help avoid confusion. 

We would like to know that the geometric description holds for < t < T, and in 
addition that properties BSO 1-8 listed below hold for all < t < T. To show this, 
we do a bootstrap argument. By IDA 1-11 and continuity of the flow u{t) in H^, 
we know that for some ti > 0, Lemma [5TT] applies giving {X(t),'y(t),r(t), z(t),b(t)) 
and e such that ORTH 1-4 hold on < t < ti with initial configuration A(0) = Aq, 
etc. Again by the continuity of u{t) (and e(t)) in we know that BSI 1-8 hold 
on < t < ti by taking ti smaller, if necessary. Now take ti to be the maximal 
time for which Lemma [5.11 applies and BSI 1-8 hold on < t < ti. (By the above 
reasoning, we must have ti > 0.) Under these hypotheses, we show that BSO 1-8 
hold on < t < ti. Since these properties are all strictly stronger than those of BSI 
1-8, we must have ti = T. 

We now outline this bootstrap argument. We have the following "bootstrap inputs" 
that we enumerate as BSI 1, etc. We assume that all the following properties hold 
for times < t < ti < T. 

BSI 1. Proximity of r(t) to 1, or localization of the singular circle 

\r{t) - 1| < (a*)^ 

and proximity of z{t) to 

\z{t)\ < {a*) 



BSI 2. Smallness of b{t), or closeness of Qb to Q on the singular circle 

< b{t) < (a*)i 
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BSI 3. Control of e{t) error by radiative asymptotic r^^t) 



^{'t)= / \^{r,z)^{r,z,t)\'^iix{t),r{t)drdz+ / \e{f, z,t)\'^e drdz <Tt^y 
BSI 4. Control of the scaling parameter X(t) by the radiative asymptotic Tfj^^t) 

x'(t)\Eo\<rl^,y 

BSI 5. Control of the r-localized momentum by the radiative asymptotic 



Im / Vi/j ■ Vu{t) u{t) dxdydz 



< r2 



X{t) 



BSI 6. Control of the scaling parameter X{t) by b{t) 

0<A(t)<exp(-exp^ 

BSI 7. Global bound on u{t) 



BSI 8. Smallness of u{t) outside the singular circle 

||Xi(rMt)||^i <(a*)i 

Assuming that properties BSI 1-8 hold for < t < ti < T, we prove that 
following "bootstrap outputs," labeled BSO 1-9 hold. 

BSO 1. Proximity of r{t) to 1, or localization of the singular circle 

|r(t) - 1| < 

and proximity of z{t) to 

\z{t)\<{a*)l 

BSO 2. Smallness of b{t), or closeness of Qi, to Q on the singular circle 

< b{t) < (a*)-^. 
BSO 3. Control of e{t) error by radiative asymptotic Th(j,) 



£{t)= / |V(,-,5)e(f,z,i)|VA(tMt)C?rdz+ / |e(f, i)| V'-^l dfdz < tJ.. 
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BSO 4. Control of the scaling parameter \(t) by the radiative asymptotic Th(^t) 

X'{t)\Eo\<Tt^,y 

BSO 5. Control of the r-localized momentum by the radiative asymptotic 



m 



Im / V^Jj ■ Vu{t) u{t) dxdydz 



BSO 6. Control of the scaling parameter \{t) by h{t) (which will imply an upper 
bound on the blow-up rate) 



< A(t) < exp ( — exp 



hh{t)) • 



BSO 7. Global bound on u{t) 

||S(t)|U. <(a*)i 
BSO 8. if^ smallness of u{t) outside the singular circle 

IIXi(r)^(t)||^. <(«*)t. 

The bootstrap argument proceeds in the following steps. The nontrivial steps are 
detailed in the remaining sections of the paper. 

Step 1. Relative sizes of the parameters A(t), T^it), a*. Using BSI 2, BSI 6, 
and ZP 2, we have 



(5.1) x{t) < r 



10 



Step 2. Application of mass conservation. Using BSI 1, BSI 3, BSI 6, and 

conservation, we obtain BSO 2 and BSO 7. In other words, mass conservation 
reinforces the smallness of b{t) and also the smallness of the norm of u. This is 
carried out in ^ 

Step 3. e interaction energy is dominated by e kinetic energy. That is, 
the e energy behaves as if it were critical and subground state. We obtain by 
splitting the term in the energy of e into inner and outer radii, using the axial 
Gagliardo-Nirenberg for outer radii and the usual 3d Gagliardo-Nirenberg for inner 
radii 

(5.2) J \e\'^fi(r)dfdz<5{a*) J | V(r,5)e|^/i(f) drdz. 

This states that in the e-energy, the interaction energy term is suitably dominated by 
the kinetic energy term. 



16 



JUSTIN HOLMER AND SVETLANA ROUDENKO 



A useful statement that comes out of this computation is 
(5.3) / \e\^drdz <\\u\\l2S{t). 

JRKlO/b 

The proof of this statement only uses the 2d Gagliardo-Nirenberg inequality (since 
we have the localization R < 10/6) and does not use the H^^"^ assumption. This is 
carried out in §H1 

Step 4. Energy conservation of u recast as an e statement. Using BSl 3, 
(15. ip . BSl 1, BSl 6, BSl 7, BSl 4, BSl 8, energy conservation, and properties of Qb, 
we obtain 



2(ei,S) + 2(e2,e)- J |V(,,5)e|V^) df-dz 
+3 / Qhldfdz+ [ Qhjdfdz <tI-^'^ + Sia*)£it). 

J\R\<lO/b J\R\<lO/b 



It results from plugging the representation of u in terms of Qb{t) and e into the energy 
conservation equation for u. The result is basically 

• the energy of Qb{t), which is small and shows up on the right side as the ^b{t)^^ 
term, 

• the energy of e, which shows up as the — J | V(f,i)epyu(f) drdz term on the left 
(from Step 3, the interaction component of the energy is small and is put on 
the right), 

• cross terms resulting from lu]"^ which are linear, quadratic, and cubic in e. The 
linear terms are kept on the left as 2(ei, S) + 2(e2, &)■ The quadratic terms 
are kept on the left as well, while the cubic term is estimated away. 

This is carried out in ^ 

Step 5. Momentum control assumption (BSl 5) recast as an e statement. 
Using BSl 5, (15. ip . properties of Qb, BSl 1, BSl 2, we obtain 

(5.5) |(e2,V(,,,)S)|<5(«*)^(t)i/^ + r^ 

The term that we keep on the left side comes from the cross term. This is carried out 
in m 



Step 6. Application of the orthogonality conditions. In this step, the orthog- 
onality assumptions are used to deduce "laws" for time evolution of the parameters 
X{t), b{t), 7(t). Using the orthogonality conditions ORTH 1-4, BSl 3, and (15. ip . we 
obtain the following estimates on the modulation parameters: 



(5.6) 



+ \bs\ < cS{t) + T 



l-Cr) 
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(5.7) 



+ 



{Ts, Zs) 



1/2 I pl"C7? 



A 

where, we recall, = —A + 1 — 3Q^. This is carried out in ^ ITTl 

Step 7. Deduction of BSO 4 from (EH). Using ([5lD, BSI 4, BSI 6, properties of 
Fb, IDA 3, we obtain 

^{X'e'-/') < 0, 
as 

which, upon integrating in time, gives BSO 4. This is carried out in ^ IT2l 

Step 8. Deduction of a local virial inequality. Using (15. 4p . (15. 6p . (15. 7p . the 

coercivity property, and the spectral property for d = 2, the orthogonality conditions, 
we obtain a "local virial inequality" 

(5.8) bs>6oS{t)-rl-f\ 
This is carried out in §131 

Step 9. Lower bound on b{s) ( =^ upper bound on blow-up rate) Using 
(ESD, (ES]), BSI 4, BSI 6, BSI 2, and IDA 1, IDA 4, we obtain 

(5.9) b{s) > 



4 log s 



(5.10) A(s) < VAoe^3isi7 

which, together imply BSO 6: 

(Km + ( 1510|) =^ BSO 6 ^ ^ < 6(t)log|logA(t)|. 

5 

This will later be used to give an upper bound on the blow-up rate. (l5.9p -f l5rT0|) are 
consequences of a careful integration of (15. 8p and an application of the law for the 
scaling parameter (15. 6p . This is carried out in 

Step 10. Control on the radius of concentration. Using (15. 7p . (I5.10p . IDA 2, 
IDA 4, IDA 1, we obtain BSO 1. This is carried out in ^ IT5l 

Step 11. Momentum conservation implies momentum control estimate 
(BSO 5). Using BSI 3, BSI 7, proof of BSO 6, BSI 4, BSI 6, IDA 3, we obtain 

BSO 5. This is carried out in 3161 



Step 12. Refined virial inequality in the radiative regime. Here we prove a 
refinement of (15. 8p in the radiative regime. Let 03 (-R) = 1 for i? < 1 and 03 (-R) = 
for ^ > 2 be a radial cutoff. With (b = 03(^/A)Cfe, A = e^"/'' for some small constant 
< a <^ 1, we define i = e — (. In this step we establish 

(5.11) dsMs) > 6iS{t) + cTb-^ [ \e\^dfdz, 

^1 J A<R<2A 



18 JUSTIN HOLMER AND SVETLANA ROUDENKO 

where 

Ms) ^ IwRQbWa^ + ^^^/ ■ VC C + (e2,ACre) - (ei,ACim) 

and 

^{t) = / \^ if,z)e{r, 2,1)1"^ nxit),r{t) dfdz + / \e{r, z,t)\'^e~^^^ drdz. 
This is carried out in ^ IT71 

Step 13. dispersion at infinity in space. Let (l)i{R) be a (nonstrictly) increas- 
ing radial cutoff to large radii. Specifically, we require that (p4,{R) = for < i? < | 
and MR) = 1 for ^ > 3, with i < ^R) < i for 1 < ^ < 2. We next prove, via a 
flux type computation, an estimate giving control on the term X4<|^|<2yi drdz: 

This is carried out in ^ IT8l 

Step 14. Lyapunov functional in H^. By combining (15.111) and (15.121) . we define 
a Lyapunov functional and show 

(5.13) dsJ <-cb(Tb + £it)- 



A<\R\<2A 

where JT" (defined later) can be shown to satisfy 

(5.14) \J-dob^\<Sib^ 

for some universal constant < (5i ^ 1, and a more refined control 

^l-Ca I 1 
b + C 
t1— Ca I /1 2 



(5.15) j(s)-hm) 
with /2 given by 



and 



/i(6) = ^||i?QJi2 + ^1^1 (r, 5) ■ VC Cdfd~z. 



< +00, and hence, (I5.15P refines ( 15. lip . This is carried out in ^ |T9 



Here, < ^ 

Step 15. Deduction of control on £{t) and upper bound on b ( =^ lower 
bound on blow-up rate). By integrating ( I5.13P and applying ( 15.140 . we prove 

An 

(5.16) b{s) < 

3 logs 
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and 

(5.17) [ (Tb^^)+S{a))da <ca*. 

Jso ^ 

Using f l5.13p and flS.lSp . we prove BSO 3, the key dispersive control on the remainder 
term e. This is carried out in §201 

Step 16. H^^"^ interior smallness. Using a local smoothing estimate and f l5.17p . 
we prove BSO 8. This is carried out in §211 

This concludes the outline of the bootstrap argument, and we now know that BSO 
1-8 hold for all times < t < T. The remainder of the proof of Theorem II .11 is carried 
out in §22H2l 

6. The equation for e 

Recall (13. ip and write 
(6.1) 

1 ^ fr- r(t) z - z{t)\ .^(,) J_ f r -r(t) z - z{t) \ • 



uit, r, .) = ^Q,(, j e-i + tj e 

In the remainder of this section, Q will always mean Qb{t)- Recall 

r — r(s) _ z — z(s) 

T = 7 = — — 

A(.) ' \{s) ' 

and thus, 

r% ~ „ _ Zg -A^ 

9.r = -y-ry, dgz = ---z-. 
Direct computation when substituting (16. ip into (13. ip gives 



\ X X J ^ \ X X 

.„ .A^ . / Tg -A^A . I Zg ^Xg . 
+ idgt - i—e + I I -— - I f^r-e + « I -y - z— ] 4-e - 7^6 

-i{^ + h^ KQ - ^(r„ Zg) ■ V^r,z)Q + ibgdbQ - ^gQ + {ihKQ - Q) 

+ idgd - i ^y + Ae- ^(r^, Zg) ■ V(f,z)(^ - 7se + {ibAe - e). 



where 7(5) = 7(3) — s. 
Also, 

and 



A^e '^A^r,z)U = A^f,z)Q + A(f,z)e 



X^e '■'^-drU = -dfQ + —dfe. 
r n n 
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Finally the nonlinear term: 

AV*^|m|\ = |Q + ep(Q + e) 

= \Q\^ Q + 2{Re Qe)Q + \Q\\ + Q\e\^ + 2{ReQe)e+\e\h 
^ " ^ ' 

linear quadratic cubic 

= \Q\^Q + [(2S2 + |Q|2)e, + 2See2] + ^[(262 + \Q\^)e2 + 2Seei] + i?i(e) + iR^ie), 

where the quadratic and cubic terms are put into -R(e). 

Adding up all of the above in fl3.ip . and taking the real part, we get: 

(6.2) 6,9,9 + 9,e2 + M+(e) + 6Ae2= (^y + 6^ A0 + ^(r„ z,) ■ V(,,5)0 - 7sS 

+ (^y + bj Ae2 + ^(r^, z^) ■ V(f,i)e2 - 7^61 

where 

(6.3) M+(e) = ei - A(,,,-)ei - -dfe^ - \ + 1 j \Q\'e, - 2See2. 
Taking the imaginary part, we get 

(6.4) 6,9feS + 9,ei-M„(e) + 6Aei = (^y + 6^ AS + i(r„ 2,) ■ V(,,i)S + 7,6 

+ (^y + bj Aei + ^(r^, z,) ■ V^f^^)^! + ls^2 
+ Im§-i?2(e), 

where 

(6.5) M_(e) = £2 - A(,,,)e2 - -d.t^ - f + 1 j \Q\\2 - 2S0ei. 

7. Bootstrap Step 2. Application of mass conservation 

In this section we prove BSO 2 and BSO 7 are consequences of the conservation 
of u{t) and several other BSI's. The assumed smallness of 6(t), e{t), and A(t), initial 
smallness of m (at t = 0) combined with norm conservation for u{t) reinforces the 
smallness of u{t) in for < t < ti and smallness of h{t) for < t < ti. 

Recall (16. ip and denote 
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Substitute (16. ip into the mass conservation law 

J \u{r, z,t)\'^rdrdz = \\uo\\'j^2 

and rescale to obtain 

Ikollia = j \Qb{t)if,z)\'^IJ^if)drdz 

(7.1) +2 Re J Qb(t){r,z)e{r,z) ^{r)dfdz 

+ / \u{r, z,t)\'^rdrdz. 



In the first term, write out /i(f) = Xr + {r{t) — 1) + 1, which splits the integral into 
three pieces. The third piece we keep; the first two we estimate using 

[BSI 6 : A < exp(-exp7r/106(t)) and BSI 2 : b{t) < (a*)^/^] ^ A < {a*Y^^ 

and BSI 1: \r{t) - 1| < (a*)^/^. The second term in (^A^) we estimate using Cauchy- 
Schwarz and the assumed control on e{t) given by BSI 3, ZP 2, and BSI 2: S{t) < 
rb(t) — {<^*y^'^- Collecting the results, we now have 



j \Qbit)\'^dfdz + \\u\\l2 = \\uo\\l2 + 0{{a* 



Finally, we use QP3: 

dib"^) = ^IIQJ^ , do = d(0) > 0, d differentiate, 

which gives 

\m\l. =\\Q\\h +dob' + 0{b' sup \d'{a)\). 

0<a<b2 

Since b^ < (a*)^/^, we have 

xyz rz 

and conclude by IDA 9 to get 

dob{tf + \\umh^^<c{ar^', 

which implies BSO 2 and BSO 7. 

8. Bootstrap Step 3. Interaction energy < kinetic energy for e 

Here we deduce (15. 2p (or (3.94) in [33]), which controls |e|^ by |Vep, and is useful 
later on. 

By change of variables, we have 

\e\'^fi{f)dfdz = X'^{t) / \u\'^ r drdz + X'^{t) / \ufrdrdz. 

J\r\<h J\r\H 
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For the first term, we use the (r-locahzed) 3d Gaghardo-Nirenberg estimate. For 
the second term, we use the axially symmetric exterior Gaghardo-Nirenberg estimate 
Lemma 12.11 Together, they give the bound 

<X^{t)\\u\\\ W^uWh +X\t)\\u\\% WuWh . 

II II IIL,^, "V,..|.|>1}" 

{xyz, |r|<|} 

We then apply BSI 8 (smallness of 11-1111^1/2 for |r| < |) and BSO 7 (smallness of \\'u\\l2 
globally in r) to get the estimate 

<{ar/'X\t)\\V^yMli^^- 

For axially symmetric functions, we have I^^Mp = l^ajMp + I^^Mp, so we may replace 
"^xyz by Vrz, then rescale to obtain 

< {a*y/^ j |V(?,5)e| dfdz 

J f,z 

and this is (15. 2p . 



9. Bootstrap Step 4. Energy conservation 

In this step, we prove (15.41) as a consequence of various bootstrap assumptions and 
the conservation of energy. 

Plug (16. ip into the energy conservation identity 



2X^Eq = X^ I \V (^r,z)u\'^rdrdz - ^A^ / \uWdrdz. 



The result for the first term on the right side is 

X^ j \V {r,z)u\'^rdrdz = J \\/{f,z)Qb{t)\'^l^if) dfdz 



+ 2Rej V {f,z)Qb{t) {f,z)(^ lJ^'{r)dfdz 
+ J \^{f,z)'^\'^ lJ'ir)dfdz 



1.1 + 1.2 + 1.3. 
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For the second term, we obtain 

— 2 J {ReQb(j,)eyii{f) dfdz — J \Qh(j^^\'^\e\'^ii{f) dfdz •(—quadratic 



— 2 j \e\ Re((56(t)e)/x(r) drdz •<— cubic 

— -y |e|V(^) dfdz 
= II.l + II.2 + II.3 + II.4 + II.5 + II.6. 

In all of these terms except 1.3 and II. 6, we can write out /i(f) = Af + r(t) and 
"discard" the Af term by estimation, since A(t) < ^Ifty We should thus hereafter in 
this section replace //(f) with r{t) for all terms but 1.3 and II. 6. We have the Qb{t) 
energy terms: 

We also have the linear in e terms that we combine, and then substitute the equation 
for Qb(t) to get 

r(i)-^(I.2 + II.2) = -2Rey" {A^f,z)Qbit) + \Qb{t)?Qb{t))edfd~z 

= — 2 Re j Qb{t)^dfdz — 26 Im j KQh{t) e dfdz + 2Re J ^h{t)^dfdz. 
The middle term is zero by ORTH 4, and the first term can be rewritten to obtain 
r{t)-\l.2 + II.2) = -2(E, ei) - 2(0, €2) + 2Re J ^b(t) edfdz. 

Next, for the quadratic in e terms in II, replace Qb{t) by Q and use the proximity 
estimates for Qb(t) to Q to control the error 

r(t)'^(II.3 + II.4) = -3 y Q^el dfdz - J Qhjdfdz + errors. 

For the cubic term in II, we estimate as (using that 1/6 <^ 1/A) 

\dfdz 

\ 1/2 / ^ \ 1/2 

^dfdzj i / \e\^\Q\^dfdz] . 



^(i)-^|II.5| <2 J |e|3|4(,)| 



<2( / \e 

'r<i/ioa 
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By 2d Gagliardo-Nirenberg applied to the first term (where is a smooth cutoff to 
R < 1/lOA) 



< 2 ( / \e\'^dfdz 
'r<i/iox 



1/2 



\V^f,z)[<f^{RX)e{R)]\'drdz 



1/2 



\e\'^\Q\'^dfd2 



1/2 



Note that if i? < 1/10 A, then /i(f) ~ 1. The first term above is controlled by some 
a* power (using the R restriction to reinsert a /i(f) factor) by rescaling back to u by 
BSI 7. The second term is controlled as 

A / \e\^drdz+ / \Ve\^dfdz\ 

^ 1/2 

< A||-it||j;^2 + I / \Ve\'^fi{r) dr dz 

For the third term, we use |Ofc{t)l ^ e~^/^. These considerations give 
The quartic term in e in II is controlled by (15. 2p . We next note that 



BSI 1 : \r{t) - 1| < (« 



*\l/2 



r(t) 



1. 



Collecting all of the above estimates and manipulations, we obtain (15.41) . 



10. Bootstrap Step 5. Momentum control assumption recast as an e 

STATEMENT 



In this section we prove (15.51) . Recall BSI 5: 



A(t) 



< r2 



Im j ViIj ■ yu{t) u{t) dxdydz 
A basic calculus fact states that for axially symmetric functions / and g, 

'^{x,y,z)f ■ '^{x,y,z)9 = {r,z)f ' '^{r,z)9- 

Plug in (16. ip into the left side of BSI 5 and change variables (r, z) to (f, z) to get 
X(t) Im Vij- Vu{t) u{t) dxdydz 



Im 



{V(r,z)^p)ir,z) ■ [V(f,z)iQb(t) + e)) ir,z) {Qb(t) + e){r, z) ij,{r)dr dz 
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then using that ipi^r) = (1, 0) on the support of Qb{t), we continue as 



Im j dfQb(t) Qb{t) fJ'{r) dfdz 
+ 2 Im y dfQb{t) ^l^'ij') dfdz 



+ Im y dft e ^{f) drdz 
= I + II + III. 

In term I, we expand out /i(f) = A(t)f + r{t) which gives two terms: the first of these 
we estimate out and use A(t) < rb(t), and the second of these is by QP 3. In term 
II, we expand out fi{f) to get 

(10.1) II = 2A(t)Im j dfQb(t)^rdfdz + 2r{t){dfQ,ei)-2r{t){dfJ:,e2). 

The first of these is estimated away using A(t) < T^j^^-j, the third term we keep, and 
for the second we note 

||e^/2(9^e||i2 ^0 as 6^0 

(r,z) 

by QP 1, and thus, we can estimate {dfQ, ei) by Cauchy-Schwarz. For term III, we 
first estimate by Cauchy-Schwarz to get 

\ 1/2 / r \ 
ep/i(f) dfdz 



liiil<(/i*.lV««i) (/ 



For the second factor, we convert e back in terms of u: 

e(f, z) = \{t)u{X{t)f + r{t), \{t)z + z{t),t) 

to obtain 

j |ep/i(f) dffiz = ||£t||^2, 
which is < a^/^ by BSO 2. Combining all of the above information, we get 
(10.2) r(t)(9,S,e2) < 5{a*)S{tf/^ + TI^^. 

We finish by using that 



BSI 1 : Wit) - 11 < (a*)^/2 ^ J_ < 2 



and then (15.51) follows. 
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11. Bootstrap Step 6. Application of orthogonality conditions 

11.1. Computation of A^/A + b. Here we explain how to obtain (15.61) . 

Multiply the equation for ei (16. 4p by R^T, and the equation for £2 (16. 2 p by R^Q and 
add. We study the resulting terms one by one. 
Term 1. 



:ii.ii 



Using the QP properties: 



-)■ as 6 -)■ 



and 



C3 



C2 



as 6 



to show that (fTTTD ~ 6, Re(iii?2g^ ^2g) ^ ^^^^ specifically, we can show that 
(HHD is 

< 1 6, 1 (5(6), where 5(6) ^ as 6 ^ 0. 
Term 2. Using ORTH 1 condition, we have 
(9,ei,i?2S) + (9,e2,i?'e) 



:il.2) 



-6 



(ei,/?2S) + (e2,/?'e) 



(ei,/?25feS) + (e2,^'<9fee) 



(ei,^2^bS) + (e2,i?'9fee) 



which we then estimate as 
< I6J 



1/2 



\t\^e-^drdz < \b,\£{ty/\ 



Term 3. 



[11.3) 



-M_(e) + 6Aei, ^^s) + (M+(e) + 6Ae2, ^^6) 
= -(e2,i?'S) + (ei,i?2e) 

+ (A(,,,)e2,i?'S) - (A(,,,-)ei,i?2e) 

+ -(9,e2,^'S)--(9,ei,i?2e) 



+ 



\Qb 



+ 1 |g,|S,i?^s 



2S2 



+ 1 |g,|^ei,i?^e 



+ (2Seei,^^S) - (2see2,^^e) 

+ 6[(Aei,i?2s) + (Ae2,i?20)] 
I + II + III + IV + V + VI. 
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By integration by parts, 

I + II = -(e2, + (ei, R^Q) + (e2, A^^M^^^)) " (^i' ^if,z){R^O)). 

Computation gives, for any function /, that A(r = 4A/ + R'^A(^f^z)f- Thus, 

using ORTH 4, 

I + II = Im(e,i?2(-4 + Agb)). 
Substituting the equation for Qb gives 

I + II = Im(e, R\-ihAQb - Qb\Qb? - ^b))- 
By examining IV+V and rearranging terms, 

IV + V = Im(e,i?24|4|2)^ 

and thus, 

I + II + IV + V = Im(e, R\-ibAQb - ^t)). 
In fact, adding VI cancels the middle term: 

I + II + IV + V + VI = - Im(e, R^-^b))- 

Properties (I3.5P and ZP 2 imply 

(11.4) \\Piy)¥'\y)U^<Cp,yf"^'\ 
thus, we get 

\{e,R'^b)\ < rl^'"^"^ ( [ |6|V^«V^' = rl'-'-'^'^sity/'. 



\jR<W/b / 

For III, use A < P^" and yU ~ 1 and Cauchy-Schwarz to obtain 

1/2 

|III| < Pb° ( / \e\'^e-^ dfdz 

\jR<W/b 

Together we have 

1/2 



mM\<rl^'~'''^ ( [ \e\'e-^drdz) 

\jR<10/b J 



Term 4- 

(11.5) (^y + [(AE, r^j:) + (Ae, ^^e)] 

By integration by parts, for any function /, we have (A/, R^f) = — J R^P-, and thus, 
the above expression becomes 
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Term 5. 



(11.6) (^^(r„ Zs) ■ V(r,z)^, R'J:j + (^-(r„ z^) ■ V(,,,)e, R'Q 

Each of these two terms, it turns out, is zero. For example. 



2 

= 0, 

since it is the integral of an odd function (recall f goes from —oo to +oo) 
Term 6. 

(11.7) 7.[(0,i?'S)-(S,i?2e)] = o. 

Term 1. 



(11.8) + [(Ae2,i?^e) + (Aei,i?2E)] < 

by Cauchy-Schwarz. 
Term 8. 

(11.9) ■ V(,,,)ei, Ffj:^ + ■ V(,,,)62, ^^e 

we estimate by 

I (r„ 2;.) I ^^^^1/2^ 



A 

similar to Term 7. 
Term 9. 

(11.10) 7.[(e2, i?'S) - (ei, R^Q)] < \is\S{tf'^. 
Term 10. 

(11.11) (Im^,^2j^) - (Re^',^20) _ lm(^fe, J?2g^). 

It turns out that this term is merely zero (< S{t)), which is shown by substituting 
the equation for in terms of Qb into fill. lip and applying the property (13. 9p . i.e., 

Im j (f, z) ■ V(f,z)Qb Qbdfdz + ^WRQbWh = 

and 

-{drQb,R'Qb)<rl'>. 
A* 

Term 11. 

(11.12) (i?i(e),i?2e)_(i?2(e),i?2s) 
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Recall that -Ri(e) and -R2(e) consist of quadratic and cubic terms. A typical quadratic 
term has the form 



R<10/b 



For the cubic terms, we use Cauchy-Schwarz, (I5.3p and properties of Qfc- A typical 
term is 



lei-^lQbldrdz < ( / \e\'^drdz 

'R<lO/b 
'R<W/b 

<S{t). 



nj 



2|,n |2 



e\-\Qb 



1/2 



1/2 



e e 



1/2 



Collecting the above estimates on terms flll.ip -( ril.l2p . keeping only flll.Sp . we obtain 
(11.13) 



+ b 



6{b)\h\+£{ty/' 



X 



+ + |7.l ) + ^(t)^/^rf -^"^ + s{t) . 



11.2. Computation of {rs,Zs)/X. Multiply the equation for ei (16. 4p by (f,5)S and 
the equation for €2 (16. 2 p by {r,z)Q and add. Note we will have now a vectorial 
equation and again study each term separately. 
Term 1. 



:ii.i4) 



hs [(a^S, (f , 5)S) + (9^6, (f , ~z)Q)] = hs Re(a,Qb, (f, z)Q,). 



Recall that a^Qb ~ \R^Q, and e^^/^^g^ _ 



— 7- as 6 —T- 0, hence, we have 



bs ReidbQb, (r, z)Qk) ^ Re(- R^Q, (f, z)Q) = 0, 



or similar to Term 1 in (lll.ip . ~ \hs\ 6{b) with 5(6) — ?■ as 6 — t- 0. 
Term ^. Using ORTH 2 condition, we have 



:ii.i5) 



(9,ei,(f,5)S) + (9,62, (f, 5)6) 

= 9, [(ei, (f, 5)S) + (e2, (r, 5)0)] - 6, [(d, (f , 5)5^2) + {e^, (f, 5)9,6)] 
= -bs Re(e, (f, z)dbQb) 
< \bs\S{b), 



since again (9{,Qfe ^ jR Q and Re(e, (r, z)dbQh) ~ 0. 
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Term 3. 



(11.16) 
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-M_(e) + bAei, (f, 5)S) + (M+(e) + bAe2, (f, z)e) 

= - Im(e, (f, 5)(5b) + Im(Ae, (f, 5)(5b) + - Im((9,,e, (f, 5)0,,) 

+ Im(e, (f, z)|Q6rQb) + Im(e, (f, 5) z6AQb) 

= Im(e, (f , 5) -Qfc + AQb + IQbPQfc + ^^AQ^ 

A 



+ 2 Im(e, V(f,2)Qb) + - lm{dre, (r, z)^^) 



-(e, (f , 5)^6) + 2(e, V(j=,s)Qfc) + -(9^e, (f , z)Qb) 



= Im 

= 1 + 11 + 111. 

For term I, we use the estimate (111.41) for to get 

|(6,(f,z)vi/,)|<rf-^^^f(t)V2. 

The term II (for example, the first coordinate of the vector with drQb) "we write as 

2[(e2, drT.) - (ei, 9,9)] + (ea, f AS) - (ei, f A0) 

and note that (9,6 ^ as 6 — )■ 0, thus, |(ei, 9,0)1 < (5(6)£^(t)^/^, and the term 
|(e2,9,S)| < (5(a*)f (t)^/2 ^ ^Yiich we estimated as in ffTTOD . The last two terms 
are estimated out by Cauchy-Schwarz and properties QP. 

For III, use A < and /i ~ 1 and Cauchy-Schwarz to obtain |III| < V^£{t)^/'^. 
Together we have 

KUni < rl^'^^^^f (t)i/2 + (5(6) + 5{a*))S{tY'^ + + v'^sitf'^ 



Term 4- 
(11.17) 

'A. 



A 

Term 5. 



(11.18) 



+ 6 [(AS,(f,5)S) + (Ae,(f,5)e)] 



A. 



M Re(Agfe, (f , 5)g6) = 0. 



^(r^, Zs) ■ V(f=,5)S, (r, 2;)S^ + ^('^s, ^^) " V(f,5)9, (r, 2;)e 

('^s, ^s) 



Re 



A 



&llL2{f,i)- 



Term (J. 
(11.19) 



7s[(e,(f,5)s)-(s,(f,5)0)] = o. 
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Term 7. 
(11.20) 



+ [(Ae2,(r,z)e) + (Aei,(f,z)S)] < 



A 



applying ORTH 2 and by Cauchy-Schwarz. 

Term 8. 

(11.21) 

^'-^ ■ V(,,,)6i, (r~, ~z)T>i + (^^^ ■ V(,,,)62, (f, 5)e) < 



similar to Term 7. 
Term 9. 

(11.22) 

Term 10. 

(11.23) 



7.[(e2,(r,5)S)-(ei,(f,5)e)]< |7.|^(t)^/'. 



(Im ^, (f, 2)S) - (Re ^, (f , 5)9) = Im(^'b, (f=, ~z)Q^). 



Substituting (13. 7p for \E';, in terms of Q;,, one can see that all terms are zero (integration 
by parts or by degeneracy of the momentum QP 3, first property in (13. 9p ) except for 
^ Im(9rOf)) (^5^)0?)) which is bounded by V]^ since A < V]^ (and localization of Qh 
implies /i ~ 1 and boundedness of the inner product). 
Term 11. 



(11.24) 



(i?i(e),(f,5)e)-(i?2(e),(r,5)S) <£(t), 



in the same fashion as Term 11 m. (I11.12p . 

Collecting the above estimates on terms f lll.14p -f lll.24p . keeping only f lll.lSp . we 
obtain 

:ii.25) 



K!ilM||g,|||,<5(5)|6,|+^(t)V2 



A 



11.3. Computation of hs- Multiply the equation for ei (16. 4p by — A0 and the equa- 
tion for 62 fl6.2p by AS and add. 
Term 1. 

(11.26) hs [(5^6, AS) - (SfcS, A0)] = h, Im(9,Qb, KQi,). 

Recalling that d^Qi, ~ \R^Q and AQb = + RdnQb, we estimate flll.26p as 



1,, ^ 
4' 



6,im(a,g,, Ag,) ^ ^||i?Q||i.(,,5) + ^ / ^^ga^g^) ^ l&.l I|i?g|li2(,,,). 
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Term 2. Using ORTH 4 condition, we have 
(9,e2,AS)-(9,ei,Ae) 

= ds [(62, AS) + (ei, AG)] + K [(ei, A9feS) - (ea, kd,Q)] 

(11.27) 

= -bs Im(e, AdbQb) 
< \bs\£itY/\ 

by Cauchy-Schwarz and using properties of Qb (e.g., see (95) in [28]). 
Term 3. 

(M+(e) + 6Ae2, AS) + (M__(e) - bAe^AQ) 

= Re(e, AQb) - Re(Ae, AQb) - - Re(a^e, AQb) 

^^^■^^^ - Re(e, IQbl^AQb + 2|g,|(SAS + GAG)) 

+ 6[(Ae2,AS)-(Aei,A0)] 
= I + II + III + IV + V. 

For term III, use A < r^°, localization of Qb, thus, ~ 1, and Cauchy-Schwarz to 
obtain |III| < r^°^(t)^/^ For term IV we use ORTH 4 to obtain 

IV = 6[((f , z) ■ Ve2, AS) - ((f, ^) ■ Vei, AG)] = Re (-2 6((f , z) ■ Ve, A4) 

For terms I, II and IV we use QP 5 scaling invariance property fl3.12p : 

I + II + IV = - Re(e, 2{Qb - ibAQb - -^b) - A^b - ibA^Qb). 

Terms with '^b we estimate by flll.4p . terms with ib we combine with the term IV 
above and estimate using BSI 2 (smallness of 6), Cauchy-Schwarz and localization 

properties of Qb, thus. Terms 3 is bounded by F^''^ '^^^ EitY/"^. 

Term 4- 

(11.29) + bj [(AG, AS) - (AS, AG)] = 0. 



Term 5. 



^(r^,2;J ■ V(f,i)G, AS j - Q(r,, ■ V(f,2)S, AG 



Im (^{rs, Zs) ■ V(^f,z}Qb, AQftj 



by the degeneracy of the momentum QP 3. 
Term 6. 

(11.31) -7s[(S,AS) + (G,AG)] = 



-7,Re(Q;„Ag,) = 0. 
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Term 7. 
(11.32) 



(^^ + [(Ae2, AS) - (Aei, AG)] = (y + ^) ^^(Ae, AQb) < 



1/2 



by Cauchy-Schwarz and closeness of Qf, to Q and properties of Q (as in Term 2 above). 
Term 8. 



X 



■ V(f,i)e2, AS 



X 



V(f:,i)ei, AO 



Im 



{rs,Zs) 
X 



V(f,i)e, AQb ) < 



7. 



:[(ei,AS) + (e2,Ae)]<|7s|^(t)'/'. 



(11.33) 

Term 9. 
(11.34) 
Term 10. 

(11.35) (Re§,AS) + (Im^',Ae) = Re(^'fe,AQb), 

which is estimated by 5(q;*), since AQtl < e^*"/'^', see Lemma 4 in [28]. 
Term 11. 

(11.36) (i?i(e), A0) - (i?2(e), AS) < £{t), 

estimating quadratic and cubic in e terms similar to Term 11 in f lll.l2p . 

Collecting the above estimates on terms flll.26p -( lll.36p . keeping only flll.26p . we 
obtain 
(11.37) 



\hs\ WRQWh < m'^' 



+ 



X 



+ l7s| + l&s 



£(t)i/2r|('-^^) + £(t). 



11.4. Computation of 7^. Multiply the equation for ei (16.41) by A^G and the equa- 
tion for £2 (O by -A^S and add. 
Term 1. 



(11.38) hs [(9,S, A^e) - (9,0, A^S)] = -hs Im(9fcQ,, A^Q,) < 5(a*) 

by the properties of (see the last estimate in Lemma 4 in ^28j). 
Term 2. Using ORTH 3 condition, we have 

(9,ei,A2e)-(9,e2,A2S) 

= ds [(ei, A^e) - (62, A^S)] - b, [{e,,A'd,e) - {e^^A'd.E 

< \hs\S{tf'^ 

by the estimate P{R)-^di,Qb{R)\ < S{t)^/^, < m < 2, from Lemma 4 in 



11.39) 
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Term 3. 

(-M_(e) + 6Aei, A^e) - (M+(e) + 6Ae2, A^E) 

+ (A(,,5)ei, A^S) + (A(,,5)e2,A2e) 
+ -(a,ei,A2E) + -(a,e2,A^e) 

(11.40) 

+ ((262 + ig.He^, A^e) + ((2S2 + ig.Hei, A^s) 

+ (2Eee2, A^E) + (2Eeei, A^e) 
+ 6[(Ae2,A2E)-(Aei,A2e)] 
= I + II + III + IV + V + VI. 

For terms III and VI we use the smallness of A and 6, localization of Qb and 
Cauchy-Schwarz to estimate them by r^°£(i)^/2. 

In terms I, II, IV and V we collect separately terms containing ei and €2- Recall 

the closeness of Qh to Q and that Q is real, thus, for example, (9^6 ~ as 6 ^ 
(recall QP 1), and hence, the terms containing drQ will be on the order of S{b) - these 
are the terms containing 62- The terms with ei produce 

(ei, -A^E + A{A''Q,) + A^^Qbl'' + 2E2) + 2EeA2e) = (ei, ^^(A^g)) + 5{b), 

where L+ = — A + 1 — SQ^, by property QP 1. 
Term 4- 

A 



:il.41) 



^+bj [(AE,A2e)-(Ae,A2E)] 
A. 



A 



+ b] (AE, (f, 5) • V(Ae)) 



< 



A +^ 



again by QP 1, closeness of Qb to Q (e.g., the terms such as drO f=:i as 6 — > 0). 
Term 5. 



;ii.42) 



A 



r., Zs) ■ V(,,,-)E, A'e ) - ( -(r„ Zs) ■ V(.-,,)e, A'E ) < 



IO-s...)| 

A 



by QP 1 similar to Term 4- 
Term 6. 

(11.43) 7s[(S, A^E) + (6, A^e)] = 7s(||Ag||i. + 5{b)), 

by integration by parts and closeness of Qb to Q. 
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Term 7. 
(11.44) 



^ + 6) [(Aei,A2e)-(A62,A^S)]< 



by Cauchy-Schwarz and properties of Qb- 
Term 8. 



11.45) 



X 



V(,,5)ei,A20^ - 



X 



■V(,,,)e2,A2s) <^^£(^)V^ 



similar to Term 7. 
Term 9. 



(11.46) 
Term 10. 
(11.47) 
Term 11. 
(11.48) 



7.[(e2,A2e) + (ei,A2S)]<|7.|£:(t)^/'. 
(Im^'^^A^e) + (Re^fc.A^S) = Re{^>b,A^Qb) < 5(«*). 
-(i?i(e),A2s)-(i?2(e),A2e) <^(t), 



in the same fashion as Term 11 iw flll.l2p . 

Collecting the above estimates on terms (lll.38p - (lll.48p . keeping only (111.431) and 
the estimate for (111.401) . we obtain 



11.49) 



7. - 



IIAQII 



(ei,L+(A^Q)) 



L2 



< £{tf" + 



X 



X 



We finish this section by observing that solving the system of equations (111.131) - 



ffTT^ - ffTT^ -f lTri^ for parameters (6„ f + 6, ^^,7^) gives (ESD and (EZD- 

12. Bootstrap Step 7. Deduction of BSO 4 from (15.61) 
From BSI 3: £{t) < ^l^t^j, we deduce from (15.61) the estimate 



\bJ < r 



1/2 

b{t)- 



Direct computation gives 
d 



ds 



(A^e^-/^) = 2A2e^^/^ 



A., 57r6, 



A 262 



12.1^ 



2AV-/M'^ + 6-6-S 
A 2o^ 
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Using that 

and (from ZP 2) 



we get that (112. ip imphes 
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< r 



1/2 

m 



262 



r 



1/2 



b(t) -ir/Ab 
- ■ l2 - ^ ' 



as 

and, integrating in s, we get 

Since ZP 2: e'^^/^^^^*)) < r6(t) < 6^3^/(4 ^W), we have 

^^^^ < X\t)\Eo\e'-"^'^ < \^mE,\e'-"^'^ < \'{0)\Eo\T;^^^' < 1 
by IDA 5, which gives BSO 4. 



13. Bootstrap Step 8. The local virial inequality 



Multiply the ei equation (16. 2 p by — A6, multiply the €2 equation by AS, and add. 
We analyze the resulting terms one at a time. 



i = 6,[(-a6S,Ae) + (a,e,AE)]. 



Since, for b small, d^Qt 
above is essentially 



by integration by parts. 
The next term is 



Q, we have d^T. ^ and dbQ ~ ^Q, and thus, the 



4 



bsiR'Q,AQ) = -b,\\RQ\\l, 



II = -(9,ei,Ae) + (a,e2,AS) 

= 9,[-(ei, AG) + (£2, AS)] + 6,[-(ei, A9,e) - (ea, A^.S)] 

which, by ORTH 4 is 

= bs[-{e,,AdbQ) - (e2, Aa,S)] < |6,|^(t)^/2. 

We continue, following the proof of [271 Lemma 8] or |28l Prop. 3]. Several terms 
come from (M_(e), A9) + (M+(e), AS) by taking adjoints using QP 5. 
Terms requiring special attention in our case are: 
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[1) The local linear and quadratic in e terms 

-2(ei, S) - 2(e2, ^) ' ^ j Q'4 - j Q 



These are related to the nonlocal quadratic-in-e term J |Vr,2ep/i(f) df dz, but 
with favorable sign, via energy conservation (15. 4p . There remain quadratic 
terms 

j fQ%Qel + j fQ^dfQel 

that are combined with j | Vf,2ep/i(f) df dz to form H{e, e), to which the spec- 
tral hypothesis can be applied. The H{e, e) term is nonlocal, and the spectral 
hypothesis can only be apphed locally, so this term is addressed with cutoffs 
as in Raphael |33] . 

2) Terms in which ^S^e and ^drQb are paired with Qi,. These are easily estimated 
due to the A factor and the localization to the singular ring. 

ei,Re(A\l/)) — (e2, Im(A\l/)). These terms are local and thus treated as in 



(3) 



14. Bootstrap Step 9. Lower bound on b(s) 
Using (15. 8p (the local virial inequality), we obtain 



(e4. 

ds 



3n 
< 1 



by ZP 2 in the last step. We now integrate this to obtain 



by the definition of sq in (13. 2p . Thus, (15. 9p follows. Now from the scaling law§ 

b 



A 



1 /9 TT 5 

- b - - g 



9 - A - 9 

A. 4 , 
^ --i > -b. 
A - 9 



'^The constant | is not very important, since e^/x ^ 1 for large x; it is chosen for convenience of 
calculations, e.g., in [33] it is chosen to be ^ . 
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Integrating this in s and inserting fl5.9p . we get 



-logA(s) > -logA(so) + / da. 

6 log a 



so 



But 



and thus, 



da ^ r log Q- - 1 s So 



so log*^ J so (logO-)^ log -5 log So' 



logA(s)>-logA(so) + ^^ ' '° 



3 Vlogs log So, 

By IDA 6: < Aq < exp(— exp 87r/96o) and the definition of So: So = e^'^^^''°, we get 
Ao = A(so) > Sq^''^'^, and thus, 

-logA(s) > -ilogA(so) + -^+ (-ilogA(so) ^° 



2 3 log s \ 2 3 log So , 

and since Sq ^ 1 (or '^1), the term in parentheses is positive and we have 

1 7TS 

(14.1) _iogA(s) > --logA(so) + 



2 ° ' ' 3 logs 

Exponentiating this equation gives fIS.lOp . Now we show how (15. 9p and (IS.lOp imply 
BSO 6. Recall ([ESD and (1530|) and that 

BSO 6 ^ - < b{t) log I log A(t) I . 
5 

Subtract log s from both sides of (114.11) to get 

— log(sA(s)) = — log s — log A(s) 

1 TC s 

> - logs - -logA(so) + 



2 ° ' ' 3 logs 

Since A(so) ^ 1, we have — |logA(so) > 0, and hence, 

7r S 71 s 

-log(sA(s)) > logs > 



3 log s ~ 3.1 log s 

Taking the log, we obtain 

TT / s \ 1 

log(- log(sA(s))) > log — + log > - log s. 

3.1 ylogs/ 2 

By dSD, 

3tt 



(14.2) log(-log(sA(s)))> g^^^^. 

Since s ^ 1 and sA(s) ^ 1, we have — logA(s) > — log(sA(s)), and BSO 6 follows. 
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15. Bootstrap Step 10. Control on the radius of concentration 
The inequality (15. 7p implies 

< 6{a*)£{tf'^ + Vl'^"^ < 1. 



Thus, 



A 

\r{s) - rol < / \rs{(y)\ da 

J so 

< I X{a)da 

J so 

/+00 
/ TT (T 

V AqC 3 log a 

< a\ 



where we applied (15.101) . It similarly follows that 

\z{s) — Zq\ < a* . 

16. Bootstrap Step 11. Momentum conservation implies BSO 5 

Recall that the function ip{x,y,z) is an axially symmetric cutoff in r only, and is 
independent of z. Also recall the convention that the default (no subscript) notation 
is for xyz space, i.e., A = + dy + d"^, V = {dx,dy,dz), and /(■■■) will mean 
/(■ ■ ■)dxdydz. Also note that (for axially symmetric u) we have |VmP = |V(r,2)Mp. 
Pair the NLS equation with |A^m + Vijj ■ Vm, integrate and take the real part to 
obtain 



Im j dtu{\A^u + V^-Vu) 

= Re y Au{\Aipu + ■ Vm) 



that we write as I = II + III. We first note that by integration by parts, we have 

1 



1 = -dt Vip ■ Vm u. 



Also, by integration by parts: 

II. 1 



We convert the second term in II. 1 to cylindrical coordinates to get 

II.l = -| / AV \u\^ + ^ I d'^ij \Vu\'^rdrdz+ ^ f dri) \\/u\^drdz. 
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Since d^i^ = 0, we calculate by switching to (r, z) coordinates and integrating by 
parts in this setting: 

II.2 = i j dlip{\dru\^ — \dzuf')rdrdz —'^ j drip\V {r,z)uf'drdz. 

Adding these, we get 



II 



-- / A'ij \U\'+ / d^.^) \drU\\ 



The two terms in III are easily manipulated (using integration by parts for the second 
one in the xyz variables): 



III = -- 



Alp \u\ 



Pulling it all together, we have 



Im / Vip ■ Vm u 



-- I aV H 



Thus, using the axial exterior Gagliardo-Nirenberg inequality (Lemma 12. ip to control 
the 4th power term, we obtain 

1 



dt Im / Wijj ■ Wu u 



< WVu 



n < c- 



'^^^^A2(t)' 

where in the last step we used the decomposition (16. ip and the estimate ||Ve||L2 < 1 
furnished by Lemma 15.11 Since we explicitly chose the origin of the rescaled time as 

So = 6^'"/^''°, we have 

■ .1 < 
da = s — So < s. 



A2(r) 



so 



We have, Vt G [0,ti), that 



:i6.1) 



m 



\m j Vi)- Vu{t) u{t) 



<\{t) 



Im / Vip ■ Vmo Mo 



c\{t)s{t). 



From f ll4.2p it follows tha 
(16.2) 



^10 



VsG[so,Si), s\{s)<T\^y 



Next, we have 



A. 
A 
A. 
A 



b- 



IF 



- A6e^"/^ 



By dESD, BSI 3, we have |A,/A+6| + |6,| < r^^'J. Also, ZP 2 implies that 1/6 ~ ln(r-(J)), 
and thus, {Girbs/b'^l < This shows that the first term is dominated by the second 



The power of T on the right side can be 5 or higher. 
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in the above display, i.e., that ds{Xe^''^'') < —^Xbe^'^^'' < 0, and thus, for all s > sq. 



we have A(s)e^^/''(") < Aoe^^/*". Using that T 
2), we obtain 



~ e' 



57r/6(t) 



1^ again a consequence of ZP 



Im / V'ip ■ Vuo Mo 



< A(t)e^"/''(*) 



Im / Vip ■ Vuq Uo 



< Ane^"/^« 



^ pSTr/feo-plO 
- ^ ^ bo 

by IDA 5. Now since ZP 2 implies ~ g-Wn/bo^ have that the above is < 1, or 



:i6.3) 



m 



Im / V^l> ■ Vmq Wo 



Inserting (116. 3p and (116. 2p into ( 116. ip . we get 



Imjvip- Vu{t) u{t) 



17. Bootstrap Step 12. Refined virial inequality in the radiative 

REGIME 

The goal of this section is to prove (15. lip . This is the "radiative virial estimate" 
which will later be combined with (I5.12p to obtain (I5.13p . which is a refinement (using 
the Lyapunov functional JT" in place of 6^) to the virial estimate (15. 8p . Once again, 
the spectral property is a key ingredient. The idea in (15. lip . (I5.12p . (I5.13P is that 
instead of modeling the solution as Qb + e, where Qh is supported inside radius -R < |, 
we model the solution as Qb + Cb + ^, where ( corrects the solution in the "radiative" 
region | < -R ^ e~'^l^ . There is still plenty of distance between the radiative region 
I < -R ^ e~'^l^ and the radius ^, which corresponds to a unit-sized distance from the 
blow-up core in the original scale. 

It amounts to taking the inner product of the e equation with suitable directions 
built on Q\, + C,. These computations are carried out in Lemma 6 of [23]. First, we 
proceed along the lines of Step 1-4 in the proof of Lemma 6 in [24] . Since the function 
/i(s) is supported in 5(0, 2A), far inside the \ width in i?, all interaction estimates 
are two dimensional - the factor /i(f) can be inserted up to possible correction terms 
with a A coefficient, which can easily be estimated using (15.11) (A < r^°). The result 
is (as at the beginning of Step 4 in the proof of Lemma 6 in ^24i|, equation (4.18)), 
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the bound 

dsfi > H{i,~e) - p^(ei,L+A2g)(gi,Ag) -CA^Eo 

+ (ei,AF,e) + (e2,AFi^) - (CrcAi^re) - (Cim,Ai^im) 

-S{a*)i£{s) + X'\Eo\)-Tl-^'^. 
By the spectral property, we have 

H{e, e) - L+A^Q)ih, AQ) > - CX^Eo - 

We also have 

/ i'2A \ 1/2 /.2A 

|(ei,AF,e) + (e2,AFi^)| <CT;/My^ |e|M <5,Tt + -J^ 
A flux-type computation (see the (4.20) and its proof in [21] ) 

-(Cre,AFre)-(Cin„AFi^) >Cr6. 

Combining these elements gives fIS.lip . 

18. Bootstrap Step 13. dispersion at infinity in space 

In this step, we prove f l5.12p . This gives us control of Jj^^pi^2A 1^1^ terms of the 
s-derivative of the 3d norm of e outside the radiative region. Since this involves 
estimating e on the whole space, the factor fi{r) is crucial, and we fully write out this 
step. With 



e 



2 



r z\ / r — r(t) z — z{t)\ 

A'a)^"^^ \X{t)A{t)' X{t)A{t)) 



we write 



dt ' ^ ' ' 



Air [ ^(1)4: ■ {dta,dt/3)\uf rdrdz + - — r-r Im [ {Vcpi ■ Vu)urdrdz, 
J AttXA J 

where a{t) = ^^t.)x{t) 1^ ^ Ait)x(t) ■ Using that ^ = compute 



... I fXs As\ f r-r(s) z - z{s) \ 1 /r, 
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We have 



2ds ^ \ A 



If h f f if z)\ 

-Im / (V04 ■ Ve)e/i(f)rffrf5 + - / — ■ V04 ( ) l^{f) df dz 

+ 6 + ^) ■ V04 j |ep/i(f) (if dz 



/ I A 



2AJ X 

Using that A^/A = —bs/b"^ and fl5.6p . we obtain that 

|ni| \Vf,,MR/A)\e\'fi{f)dfd~z. 

Using ( (I5.7P ^» |(rs, < 1), we obtain 

livi<AiiM|ir2 <r2 



xyz 



Term I is estimated via Cauchy-Schwarz and basic interpolation 

|I| < ^ y \V(t)i{R/A)\\e\^^i{f)dfd~z + Tl''^ j \Ve\'^^l{f)dfdz. 

For Term II, we use that 

^l2R.^^J^l23\>h for A< |(r~,5)| <2A 
^ V ^ / 1 otherwise 



and the fact that ^{f) ~ 1 for z)| ~ A to obtain 

b f 

II > — / lel'^dfdz. 
- 30 

finishing the proof of fl5.12p . 



y4<|(r,2)j<2A 



19. Bootstrap Step 14. Lyapunov functional in 

In this section, we exhibit the Lyapunov function J ^ b'^ and prove the upper 
bound on dsJ given in (I5.13p . This is obtained by combining (15. lip (virial identity for 
dynamics in the soliton core) and fl5.12p (dispersion relation in the radiative regime) 
and the conservation (a global quantity that links the two). 
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Let 



J{s)= y ig,i'-y g' + 2(ei,s) + 2(62,0) 



where 



m = \\\yQ,\\l^ + \\m y"(f,z)-VCC- 

The argument follows the proof of Prop. 4 in [2l]. Multiply (15. lip by (5i6/800 and 
sum with (15.121) to obtain 
(19.1) 



~ b 
> —S + — , 

- 800 800 A<R<2A 



\e\^dr dz + ^^^^b - ^X^Eq- T^^^ j \V i^f,z)(^? ^r) dr dz. 

The last term on the right side is estimated as 

Vfj |V(,,,)e|V(r)rfrd5<rf (^F^^" + j \V ^rrz)~A^ ^^{r)dr d~z 

< r;+"/^ + Vf j |V(,,,)e|V(r) dfd~z. 

Using that 

/i(s) = /i(s) + (e2,ACrc)-(ei,ACiin), 
we also rewrite the second term on the left side as 

hdsh = dsihh) - bsfi 

= dsibf,) - bsh - &s[(e2, ACrc) - (ei, ACim)] 

= ds{bh{b) - I h{v)dv + bie2,AU-b{ei,ACin 



&.[(e2,ACre) - (ei,ACim)]. 
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(19.2) 





f,4 






800 





> 



6lb_ 
800 



fi{v)dv + bie2, ACre) -bie^,ACu 

-2A 
'A ' 



1000 



The conservation of norm, written in terms of e and Qh is 



ef'^{r)drdz+ I \Qb\^ ^{r) drdz + 2B.e I eQ},^{f)dfdz 



which we rewrite as 



(19.3) 



e|V(r)rffd5 + r(t) |Q,|'- / \Q\^ + 2Re{e,Q 



\uo?-r{t) / Q2-2ARe(e,fQb). 



Then write 

(t)i{R/ A)\t\^ ^{f)dfdz 

and obtain by taking dg of (119. 3p 

1 



t\^^i{f)dfdz - / (1 - 04(i?/A))|e|V(^) drdz 



r[s] 



4(i?/A)|e|>(r) drdz 



rdz I 



+ a 

-2a 



14 



r(s) 
A 



g2 + 2(ei,S)-2(e2,e) 
(1 -04(^M))|e|V(^) ^^^f^^j 



r(s) 



Re(e, rQfc 



The last two terms are easily estimate with (15. 6p and (15. 7p . Inserting into (119. 2p . we 
obtain fl5:T3|) . 

We next show the two estimates (15.141) and (I5.15P for J . 

We first show (I5.14p . We have ||Q6||^2 — 1|<5||^2 = ((io + o(l))6^, and now we estimate 
the rest of the terms in the definition of J . By Cauchy-Schwarz, 



ei,S) + (e2,e)| < / |e|V«« < T^'^ « h 



.1/2 
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[1- MR/A))M'^Kr)dfdS < 



r-r(t)\<XA 
\^-z{t)\<XA 



The set | |r - r(t)| < \A,\z - z{t)\ < \A} has a volume ~ (XAf in M^. 

Therefore, by Holder, the above is bounded by 

1/2 



(XA) 



\r-r{t)\<\A 
\z-z{t)\<XA 



u\ dxdydz 



Then in the L^, we widen the spatial restriction to r > | and apply the axial 
Gagliardo-Nirenberg (Lemma I2.ip and the definition of A to obtain 

AA||?i||i2||VM||L2 < A (^j^ \V(^r,z)A'^l^{f)dfdi^ < h. 

We claim that |/i(&)| < h. Indeed, we first note that i6||i?(56|||2 ^ h. Also, recalling 
that C = (1>3{R/A)(, we have 



Im / Rdj^C (drdz 



< 



R\df,C\\C\drdz + j 



R\C\'drdz 



< ||i?CllL-^||VC||L2 + jWRCWU J R~ dfdz 

which establishes that |/i(&)| ^ b. It is similarly straightforward to show that 
Jo fii'v) dv < b. Finally, the local terms {e2,A(rc) and (ei,A(^iin) are estimated 
3y Cauchy-Schwarz and the ZP 1-2 properties, and shown to be ^ b. This concludes 

the proof of ( 15.14p . 

For the proof of (15.151) . see the proof of Prop. 5 in [21], which uses the conservation 
of energy. 

20. Bootstrap Step 15. Control on £{t) and upper bound on b{s) 



Next, we prove (I5.16p . Using that (15.141) implies 



J " 



^ 1, we have 



bn [do 



Thus, by (15131) (note that - J, > 0) 



exp 



g I b 



5tt do 



5tt do 



Js 



4 V J 



4 V J 



bV, 
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Integrating from Sq to s, we obtain 



5tt do \ 571 do , ^ 



4 V J{s) J V ^ V ^^^0 

Since exp {^\J^^^ > e'"/''" > Sq, we conclude that 



and f l5.16p follows. 

Next, we obtain f lSTTj) . Divide f lSlTSj) by J^/^ ^^^^ ([5:i4|) [ J h'^) to obtain 
Integrate from so to s to obtain 



/■ Vm^) + ^^(^)) rfor < j''\so) - J"%S) < h{so) < 
J so 



a 



where in the last step we applied IDA 2. 

Finally, we prove BSO 3. Let s G [so, si) (recall that we are carrying out the 
bootstrap argument on < t < ti, and si = s(ti)). If bs{s) < 0, then BSO 3 at s 
follows immediately from (15. 8 p at s (the local virial identity, bs > So£{s) — ^l^J^^)- 

If bs{s) > 0, let S2 G [so,s) be the smallest time such that for all a e (52,5), we 
have bs{a) > 0. Then b{s2) < b{s) and either S2 = Sq or 65(52) = 0. In either case, we 
claim 

(20.1) S{S2) < Tilly 

If S2 = So, then fl20J^ is just IDA 4. If 6,(52) = 0, then ( ]201|l just follows from the 
local virial identity (15. 8p at S2. From the second of the two estimates in (I5.15P and 
also (120. ip . we have 

By the first of the two estimates in (I5.15p . the fact that J' is (nonstrictly) decreasing, 
and the above estimate, we have 

f,{b{s)) + ^£{s)<J{s) + T'/' 



5/6 
b(s) 

5/6 p5/6 



< j{s2) + rf 



<Mb{s2))+Tii:^,+Tii:y 

Combining this estimate with 



df2 
< — 
db^ 



< +CX) and 6(S2) < b{s) ^ /2(6(s2)) < /2(6(s)) 
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and 

b{s2) < b{s) =^ rf,(,2) < rb(^) 

yields 

f2m) + ^£{s)<Mb{s)) + 2Tll^y 
Canceling f2{b{s)) yields BSO 3. 

21. Bootstrap Step 16. H^^'^ interior smallness. 

21.1. Preliminaries. We consider two fractional derivative operators in our analysis: 
-D^2 and D^^^. The operator D^^ acts directly on a function f{r',z') defined in 
and returns a function of (r, z) G M^: 

where J-2 is the Fourier transform on and we adopt the practice of writing the dual 
variables of {r,z) as (p, C)- The operator D^y^ acts on a function f{x',y',z') defined 
on and returns a function of (x, y, z) on 

where we adopt the practice of writing the dual variables of {x,y,z) as {C,, ■ 
Now if f{x, y, z) is an axially symmetric function on R'^, then D^y^f{x, y, z) is axially 
symmetric and can thus be viewed as a function of (r, z) G x M and, when extended 
to an even function in r, it can be viewed as a function on M^. In this case, we can 
write 

where Ti is a. zero-order Hankel transfornj^ (p. 341 of Bracewell [2]) in r and a Fourier 
transform in z: 

/' + 00 /"OO 

nf{p, 0= / JoMe^^</(r, z)r dr dz 

J z=—oo J r=0 

and the zero-order Bessel function on M is: 

(which extends Jo{oj) as an even function in u, and thus T-Lfi^p, () is even in p.)§ We 
note that 

;'+oo ;"+oo 

n''Fir,z)= / / Joirp)e'^'^FipX)pdpdC. 

Jc=-oo J p=0 



^We will not bother with constant factors (such as (27r)^^) in the definitions of the Fourier and 
Hankel transforms, as they are irrelevant in the analysis. 

^Note that if we wanted to write the domain of r-integration in the definition of H as (— oo, +oo), 
assuming /(r, z) is even in r, we would need to put |r| in the integrand in place of r. 
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If a function h{x, y, z) is axially symmetric and the support of h is contained in 
< 5 < r < 5^^^ then we will say that h is shell supported. We note that for a 
shell supported function h, \\h\\]j^^ ~ ll^ll-LSy^- Unfortunately, if h is shell-supported, 
then fractional derivatives -D^^/i are no longer shell-supported. To deal with this, and 
handle the conversion from D*^ to -0*^^, we use the following standard microlocal fact, 
which is a consequence of the pseudodifferential calculus (see Stein [35], Chapter VI, 
§2, Theorem 1 on p. 234 and §3, Theorem 2 on p. 237; see also Evans-Zworski [8]). 

Lemma 21.1 (disjoint smoothing). Suppose that dist(supp t/^i, supp ■?/^2) > and 
s, a > 0. Then for h{r, z) we have 

and for h{x, y, z) we have 

Next, we address the matter of converting from D^^ to D^yz- We need the compo- 
sition (Bracewell [1]) 

jr^T/-! = A, 

where A is the Abel transform (see p. 351 of Bracewell [2]) 
(which is the identity in Q. 

Lemma 21.2 (fractional derivative conversion). For any < s < 2, the composition 

rz xyz 

is hounded as an operator Lf.^ — )■ L,^^ . 

We remark that scaling requires that the weights in any such estimate be r" on the 
left and r~" on the right for some a G M. We only need a = | in our analysis and 
have not explored the possible validity of other values of a. 

Proof. Note that the composition under consideration is 

Note that (r + iOy^'^J^"^ = T^^dpiy^ where J" is the fractional integral operator 
/°/(p) = /^|^(p' - p)"'V(p') dp'. This, together with the fact that J^-^ is an L^^ 

unitary map, and the fact that p^^'^l-Lr~^/'^ is an L^^ — )■ L^^ unitary map, we 
reduce to proving the L^^* — )■ boundedness of 

u = d,ii'\{p'f + er''A{{p"f + e)-''\p'r"' ■ 
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We obtain Uf = dpiy^h, where 
By Fubini's theorem, we obtain 

fp" =+0O 



pp-=+oo 

Ufip,0^d, / K{p,p")f{p",C)dp", 
Jp"=p 



where 

Kin n"^ = r"'" 1 ( {P'f + e y'" ( P" , , 

^^'^^ ip'=p (p'-p)^/^(p"-P')^/nW + CV Vp" + p7 ^' 

By changing the variables twice, first p' ^ p' -\- p and then p' — (t{p" — p), we have 
K(nn")= C 1 / (a/+(l-a)p)^ + C^ y/V p" V'\ 

Thus, 

/• + 00 

[//(p, C) = -K(p, p)/(p, C) + / d,K{p, p")f{p", C)dp". 

Jp"=p 

We have K{p,p) = 2~^/^ jj^^ o-"^/^(l — a)~^/'^da < oo. For the second term, we 
change variable p" — pp to obtain 

r+oo 

(21.1) uf{pX) = -cf{pX)+ / p{dpKyp,pp)f{ppX)dp. 

Jp,=l 

But we have pdpK{p, pp) — I{pX/p)j where 
and 



Denote by Ik the result of substituting into the expression for I. Since | J2(/x, A) | < 
{p + 1)~^, we have 

P' 

If s < 2, we have 

'^^'^'"1=0 ^ ^ (p2 + P)f((ap+(l-a))2 + A2)i-f 

(1 - a)V2 



<-/' - 



V2(ap+ (1 



d(7 
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By separately considering the regions < a < i and I < a < 1, we obtain 

^1/2 da 



1 r'^' 

|/i(-",A)|<-+ / - 



i if s > i 

fj, 2 
< I i^^^^ if S = i 

if < s < i. 

Applying the norm to ( I21.ip and using Minkowski's integral inequality, we bound 

by 



IIW,C)I|l^< II/(^C)I|l^+ (^'^^max(/i-\./.-i->-trf/iJ ||/(-,C)IU- 

Following through with the norm, we complete the proof. □ 

The following were established by Strichartz |37], Keel-Tao [IT] and Kenig- Ponce- 
Vega [19j. We say that {q,p) is 2d admissible if ^ + ^ = 1 and 2 < p < oo. 

Lemma 21.3 (2d Strichartz estimates and local smoothing). Suppose that h{r,z,t) 
satisfies 

idth + Arzh = g. 
Then for 2d admissible pairs and {q2,P2), 



^IIl^iL^i ^ ll^ol 



r 2 



\{Drz) ^'"^gWilLi^ ifsuppg compact, 

where ho{r,z) = ho{r,z,0) is the initial data. In the case of the last bound supp^f 
should be contained in a fixed compact set for all t. 

We say that {q,p) is 3d admissible if | + ^ = | and 2 < p < 6. 

Lemma 21.4 (3d Strichartz estimates and local smoothing). Suppose that f{x, y, z, t) 
solves 

idtf + l^xyzf = g- 
Then for 3d admissible pairs {qi,pi) and (^2,^2); wg have 



who', p!, 



) ' 9\\l1l%^ ifsnppg compact, 



'^We thank Fabrice Planchon for pointing out that, although the endpoint case {qi,pi) — (2,6) 
is available for the local smoothing estimate, it does not follow from the Christ-Kiselev lemma. We 
will avoid the use of the endpoint in our analysis. 
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where fo{x, y, z) = f{x, y, z, 0) is the initial data. In the case of the last bound supp g 
should be contained in a fixed compact set for all t. 

The local smoothing estimates in the above lemmas are established in the homo- 
geneous case by [19] Theorem 4.1 on p. 54. The inhomogeneous version as stated 
above then follows from the Christ-Kiselev [3] lemma. 

Finally, we record some Gagliardo-Nirenberg embeddings for shell-supported func- 
tions. 

Lemma 21.5 (embeddings). Suppose that q{r,z) = q{x,y,z) is axially symmetric 
and shell- supported. Then 



(21-2) MLtLi,^<\\q\\/d^. Iivg|li41. , 



1/2 ||^^||l/2 

xyz 

(21.3) \\D%im,.f^,, <ikiil?,. . 

^ ^xyz ^xyz ^t^xyz 

Proof. First we prove (121.21) . We have the 2d Gagliardo-Nirenberg estimate 

hWlf, ^ lkllL2^l|Vr^9|li2^ ■ 

rz rz 

Since q and Vrzq are shell supported and iVr-zQ'l < iVj^j^^gl, we obtain 

\\(1\\l%^z ~ \\^xyzqtL2 ■ 

Integrating in time, we obtain (I21.2p . Next we prove (121. 3p . We begin by noting 

||nl/2| |2||2 _ |/r)l/2| |2 r)l/2|„|2\ I 
W^xyzm Whly^ — W-'^xyzm -.J^xyzm lxyz\ 

= {\(l\'^,D\q\'^)^y^ 

= {\q\',RV\q\').y, 

= (|gp,i?(RegVg))^yz, 

where R is the vector Riesz transform. Using the boundedness of the Riesz transform 
on we obtain 



xyz 

3 

xyz ' 



3/2 



< Ikllia l|Vg||L2 



iil/3ii ii2 /3 

By the 2d Gagliardo-Nirenberg estimate ||g||L6^ < \\q\\j^2 \\'^rzq\\i2 , the fact that q 
and Vrzq are shell-supported, and iV^^g] < |Va;y2g|, we obtain 



\\Dlil\q\'Uy^ < \\q\\'/2 \\V.yM%" ■ 

y J-^xyz ^xyz 

Integrating in t, we obtain (I21.3p . □ 
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21.2. Outline and notation. Effectively, we would like to restrict to outside the 
singular circle and reduce matters to the local theory, which is set in H^^^ . 
Let 



Xi{R) 



ioi R<\ 
for i? > I 



and set Xi(^5 z) = ~ + z^Y^"^). This is a cutoff to outside the singular circle, 

and we call the support of Xi the external region. Let 



UR) 



for i < < I 

for i? < I and i? > I ' 

forl<i?<| 

for < ^ and i? > i| 

— 16 — 16 



and let Xjii^,z) = — 1)^ + z'^Y^'^) for j = 1,2. Note that X2ir,z) is 1 on the 

support of Xi(l — Xi) S'lid Xsi'i^yz) is 1 on the support of X2{f^,z). We will call the 
support of X2 the tight singular periphery and the support of Xs the wide singular 
periphery. Let v = XiU, w = X2U and q = XsU. 
Our goal is to obtain the estimate BSO 8, i.e., 

(21.4) WDlilvhrLl,^ < {a*f\ 

This will be a multistep process - estimates on q (Step A) obtained from BSI 3 will 
imply estimates on w (Step B), which will imply estimates on v (Step C). 

21.3. Step A. Estimates on q. Recall q = xsm, i.e., u restricted to the wide singular 
periphery. Control on q is inherited from the local virial arguments in the previous 
sections. Specifically, mass conservation, fl5.17p . and BSI 3 imply (by rescaling) that 

(21.5) Ml-li^. < hohi^^ , II VglU.^i^^ < (a*)^° . 
We also obtain an upper bound on ti. Indeed, from fIS.lOp . 



si r+oo 



(21.6) ti = / X\s) ds < CXo / e-("/3)(s/iog.) < ^* ^ 

J so J2 

It will be understood that t is always confined to [0,ti) below (when writing L'^, 
etc.) 

21.4. Step B. Estimates on w. Recall w = X2U, i.e., u restricted to the tight 
singular periphery. We use the control on u in the larger wide singular periphery 
obtained above in Step A to provide control in a stronger norm on the smaller tight 
singular periphery. Since w is supported away from r = and away from r = 00, it 
effectively solves a 2d NLS equation. This effective 2d equation can be analyzed using 
the (locally-in-time) stronger 2d Strichartz estimates as well as the 2d Gagliardo- 
Nirenberg estimates. 
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Specifically, the goal of this subsection is to obtain, using (121 .Sp . the following 
bounds on w: 

(21.7) \my.Y'MLlLl^.<M"' 

for all 2d admissible {q,p)- 

The function Dl/yzW = Dlylx2U is still axially symmetric but unfortunately (due 
the nonlocality of Dxyz) no longer shell supported. By Lemma 121.11 we have that 

1 /2 

(1 — xs) Dxyz X2U is infinitely smoothing, and thus, 

\\{l-X3)Dl^lw\\^L-^^ < \\w\\l^l2^^. 

1/2 

Hence, to establish (121.71) . it suffices to prove, for p = Xs Dxyz w, the estimates 

(21.8) lbllL?LS,.<lbo||Li,. + («*)'/' 
and 

(21.9) ||^|U^,.^^<||^o||Li,. + («*)^/^ 

The equation (121.71) will then follow from IDA 7 and IDA 8. 

We begin with (121. 8p . Note that p is axially symmetric and shell supported. More- 
over, it solves an equation of the form 

idtp + ArzP = J^^Gj, 
j 

where the inhomogeneities can be put into the forms 

Gl = X3 Dl/yl X2 (idt + Arz)u 
G2=^3D'Jyli;2U 

G,=Vrz^3Dlilij2U. 

For each type of term Gj, the functions 1^2 and ip3 are axial cutoff functions with 
support in i < ((r — 1)^ + 2^)^/^ < | and | < ((r — 1)^ + 2;^)^/^ < |, respectively. 
For term G3, we substitute the equation {idt + Arz)u = —r^^dj-u — \u\'^u and then 
reexpress the term involving r~^drU to be of the type G2 and 6*3. At this point, the 
inhomogeneities take the form 

G2 = ^Ps D% ^2 u 

G^ = Vrz^^D\ll^2U. 

By Lemma [21.31 

(21.10) WpWlIli^ < IbolUp. + 11^111^4/3^4/3 + IIG^IUjLP. + ||(A..)-'/'G3|U?LP., 
since G3 has compact support in r and z. 
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Since Gi is shell supported, the rz norm converts to an xyz norm. The fractional 
Leibniz rule and the identity X2u\u\'^ = allow us to bound as follows: 



I r4/3,2 



(21.11) \\Gl\\.4/:i.4/3 < \\DIIIw\\l^]^4 ||g||r4r4 +||w^||l?°L4 H-D^izl^l 

Since w is shell-supported, we can apply the 2d Sobolev embedding 

II II -^t ^xyz II II -^t ^rz 

(21.12) < \\DIIM\l^lI^ 

< WxsDl/MlLrL^.. + 11(1 - X3)Dl/MLrL^..- 

We note that we have 

= (Xsr-^/^) (r V^Z^; V V- V2) (.V^Z^V^) . 
By Lemma [21.21 and the support properties of Xs, 

Wx^dHMltl^z ^ V"D%lw\\Lr^z = \\D%lw\\LrLi^^ 

^ Ml^LI^z + 11(1 - X3)Dlilw\\L^Ll^^. 

Applying Lemma 121.11 to the second of these terms, 

(21.13) WxsDli'wh^L^^^ < Ml-li,^ + Ml^li,^. 
On the other hand, by Lemma 121. H 

(21.14) 11(1 - x3)DyMLrL^.z ^ Miri^rz- 

Combining fim^ . fim^ . and (EUl, 

(21-15) IHlL^Lt^z ^ WPh^Llyz + MLrL'i^z- 

1/2 1/2 1/2 

Returning to fl21.1ip . writing DjyzW = x^DjyzW + (1 - XzjDxyzW, we have 

Pxi>IU^L4^^ < IblUfLi^, + 11(1 - xz)D%lx2u\\LrLt^.- 

Applying Lemma 121.11 to the second of these terms, we obtain 
(21-16) \\D%lw\\^,^4^_^ < WvWLiLt^.^ + M\LrLl^.- 

Combining (l2rTT|) . (|2LT5|) . (l2rT6|) . ([212]), fl2L3|) . and (l2T3|) . we obtain 

(21.17) 11^111^4/3^4/3 < («*)'(|b||L4L4^, + IblUrLp. + MWlO- 

Since G2 is shell-supported, we convert the rz integration to xyz integration and 
appeal to (121.51) to obtain 

(21.18) \\G2\\LlLl^<tM\LrLl^z+i'l^\\^^yMLlLl^^ < ("l'^'- 
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Finally, for G3 we take ip^ to be a smooth function that is 1 on the support of 1^2 
and ip3 but near r = 0. Then decompose 
(21.19) 

By Lemma [21.11 

(21.20) <tTMLrLl,A\\yq\\L^,Ll,^ 

Note that 

^^Dll'i^.Dllii^2q={i'.r-^^^^^^^^ 
and thus, 



By the fractional Leibniz rule 

(21.21) ||^4/^.?^3/^^.£^2g|U?L., < lkllL?Li,. + l|Vg|L.,|^^ < {a*Y/\ 

Now 
(21.22) 



Collecting ( l2T:i9|l . (|2L2Q]), (I2L2T]) . (l2T:22ll . 



(21.23) ll(^..)-^/'G3|L.,.p^ < (a*)^/^ 

Combine 021. lOp . 021. 17p . 021. ISp . and 021.23P to obtain (with conversion of some 
xyz-noim.^ on p and w back to rz norms) 

MlIl^. < {a*y^\\\p\\LtLi^ + IblU^LP. + Ml^l'^J + {a*f'' , 

from which it follows that 021.81) holds once 021. 9p is available. 

Next we prove 021. 9p . Note that w is axially symmetric and shell supported. More- 
over, it solves an equation of the form 

idtW + ^rzW = Gj, 

j 
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where the inhomogeneities can be put into the forms 

Gi = X2 {idt + ^rz)u 

G2 = 1p2U 

G3 = Vrz'ip2 U . 

For each type of term Gj, the function ip2 is an axial cutoff function with support in 
i < ((r — 1)^ + z^)^/^ < |. For term Gi, we substitute the equation {idt + ^rz)u = 
—r~^drU — \u\'^u and then reexpress the term involving r~^drU to be of the type G2 
and G3. At this point, the inhomogeneities take the form 





= X2 \u\^U 


G2 


= ^p2U 


G3 





By Lemma [21.3t 

(21.24) IM^LZy^ ^ IkolUp^ + II Gill ^4/3^4/3 + ||G2|UiL2^ + ||G3|UiL2^. 
For Gi, we use that q is shell-supported and apply fl21.2l) to obtain 

11^111^4/3^4/3 < |k||i4^4 <\\q\\%' \\Vq\\%l. <{ar. 
For the remaining inhomogeneities, we argue as follows: 

\\G2\\Lm^ + II^3||ljlp. < tM\LrLi^.+tT\Wq\ULi,^ < 

Collecting, we obtain f l21.9p . 

21.5. Step C. Control on v. We use the control on u in the tight singular periphery 
to obtain control on u in the external region. In this region, u solves a genuinely 3d 
NLS equation, and we thus need to work with the 3d Strichartz estimates, which are, 
locally-in-time, weaker than the 2d Strichartz estimates. Also, we must only use the 
3d Gagliardo-Nirenberg estimates. 

Specifically, we use f l21.5p and f l21.7p to prove 

(21.25) llDlilvh.L^^^^ < \\Dli%\\Li,^ + {a*fl\ 

Via IDA 8, we obtain f imj) . 
The equation for v is 

idtv + /^xyzV + \v?v = y^-^j, 

where 





= ■ip\u\'^u 


F2 


= ipu 


F3 


= u, 
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where in each term Fj, the function if) is supported in the intersection of the supports 
of Xi ci'iid 1 — Xi- Recall that X2 and Xs are 1 on this set. By Lemma [21.41 

+ II -Dt{7-^2 II L^L^ + II -^3 II L2r 2 . 

Note that for F^, we have used the local smoothing property, and the fact that F3 is 
compactly supported. 

For Fi, we use that Fi = ipwlwl"^, the fractional Leibniz rule, and estimate as 

^Dli^F,^^,^^ < Pi/>llLfL.,Jklli,^..^, < 
by 021. 15p and fl2L7l) . For F2, we estimate as 

II xyz ^Wlj^l^^y., J-ll\ xyz/ l^xyz r-^ \ J 

by (121.71) . For F3, we estimate as 

< IkllL^Ll^, + l|Vg||i2i2^^ 

< («*)^/^. 

We also estimate, using the fractional Leibniz rule, and the 3d Sobolev embedding, 

\\Dlil{\v\'^v)\\ 10/7 10/7 < llD^gwII 10/3 10/3 ||t;||^5^5 

<^ II nl/2 II II nl/2 i|2 

^ \\j-'xyz'V\\L^^ri ^10/3 \\J->xyz'^\\^5^Z0ni ■ 

f xyz j- xyz 

Collecting, we have 

ip^g^iu^Ls,. < WDiiivoUy^ + ii^y/.HLw3^io/3pi/>ii^5^3o/ii + («*)^/^ 

z y ^ xyz 

which yields (121251). 

This completes the proof of BSO 1-8. 



22. Finite time blow-up and log-log speed 

In this section, we deduce the finite-time blow-up and log-log speed. Specifically we 
deduce (122.61) . (I22.7p . and (122.81) below. The convergence (11.71) claimed in Theorem 
11.11 is a refined version of (122.61) that can be obtained by following the techniques of 
[211 Prop. 6]. 

From the above bootstrap argument, we conclude that ti = T, and hence, by (121. 6p . 
T < a* < 00. From the local theory in H^, we have that ||M(t)||j^i +00 as t T. 
By BSO 3 and the fact that ||VM(t)||L2 < A^^dlQ^lUa +£{t)), we conclude that 



(22.1) 



A(t) ^0 as t^T. 
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Also, since |As/A| < 1 on [so,si) (see BSO 2 and (I5.6p . recall Si corresponds to 
ti = T), integrating from sq to s, we obtain 

I logA(s)| < s - So + I logA(so)| . 

By (122.11) A(s) — > as s — j- si, thus, the above line implies si = +oo. 

We now claim that for t close enough to T, A satisfies the differential inequality 

(22.2) l.<_(A2iog|logA|)i<C 

for some universal constant C > 0. To establish (122. 2p . we first show that 

(22.3) ^ < 61og|logA| < C. 

To prove (122. 3p . observe that the lower bound is equivalent to BSO 6. For the upper 
bound, first note that (15.91) and (I5.16P give 

(22.4) b{s) 



log s 

Moreover, (15. 6p gives 

(22.5) - ^ ~ 6. 

A 

Integrating (122. 5p and using (122. 4p . we obtain 

-logA(s) + logA(so) ~ / h{a)da 

J so 

da 



so 



log a 



S So 



log s log So 

Recall that sq = e^'^/^'"' > 1 and < A < 1 by BSI 6. These imply logA(so) < 0, 
logA(s) < 0, So/log So > 0. Consequently, 

|logA(s)| < -^ + |logA(so)|. 
log s 

Taking the log of both sides, and taking s > so ^ 1 such that s/logs ^ \ logA(so)|, 
gives 

1 



log I log A(s) I < (log s — log log s) ~ log s 
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which gives the upper bound in fl22.3p . Now we proceed to justify fl22.2p using f l22.3p . 
Since A ^ 1, we have 

^(A^ log I log A|) = AAi(log I log A|) (2+ ^ 



dt \ I log A I (log I log A 1)2 

~ A At (log I logAI) 

= y(log|logA|). 

By fl22:3|) and f l223|) . we obtain (HT^ . Integrating (^T^ from t up to the blow-up 
time T (where A(T) = 0), we obtain 

A2(t)log|logA(t)|~(T-t). 

This implies 



(22.6) A(t) 



T-t 



iog|iog(r-t)| 



From the relationship between X{t) and t given by ( 122. 6p . the relationship between 
b{s) and s given by (122.41) . and the relationship between b and A given by (122. 3p . we 
obtain the relationship between s and t as 

(22.7) I log(r - t) < s < I log(r - t) for some ai < as 
and the relationship between b and t as 

(22.8) b ^ 



log|log(T-t)| 



The proof of the exact convergence in (122.61) and (122. Sp follows [2H Prop. 6]. This 
proves ([TTD, ([IS]) and ffH)]) . 



23. Convergence of the location of the singular circle 

In this section, we just note that parameters r(t) and z(t) converge. Indeed, by 
(15. 7p . we have jr^l/A < 1 and \zs\/X < 1. Hence, for < ti < ^2 < 7", we have 

\r{t2) - r{h)\ < J \rt\dt = - ^ dt < -. 

By the estimates obtained on A(t) above, this implies that r(t) is Cauchy as t — )■ T. 
Similarly, we obtain that z{t) converges as t — )■ T. BSO 1 (which we proved in 
Bootstrap Step 10) then implies that \r{T) — 1| < (0;*)^^^. This proves (II. 5p and 
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24. CONVERGENCE OF THE REMAINDER 



In this section, we will establish the existence of u* G L^(M^) such that u ^ u* m 
L2(M3). This will prove ([L3]) in Theorem [HI 

First, we observe that there exists u* G Lfjj^(R^\{0}) such that for each i? > 0, 



(24.1) lim / \u{t) - u*\^ dxdydz = Q . 



This is essentially a consequence of the fact that, as in Bootstrap Step 16, we can 
carry out the local well-posedness theory in H^^"^ in this external region. This local 
theory gives a continuity of the flow (in this external region) right up to time T. 

Having established that (124. ip holds, we note that it follows from (124.11) that for 
each i? > 0, 



\u*\^dxdydz < 2\\Q\\ 



'|r-r(T)|2+|2-z(T)12>ij2 

and thus, in fact u* G L^(]R^). We next upgrade (124. ip and show that in fact 



(24.2) u{t) u* in L'{W) 

(not just in Lj^Qj,(]R'^\{0})). The following estimate is, in the 3d straightforward 
consequence of the Hardy inequality. But in 2d, it is a little more subtle and requires 
a logarithmic correction: 

Lemma 24.1. Suppose that (j){x,y) is a function on M^. Then for any D ^ 1, we 
have (with r = y^x^ + y"^), 

\\<i'\\h(r<D) < ^'(Il0llii,(.<l) + (log/^)||V0||i.^(,<^)) ■ 

Proof. We follow [2^ Apx. C] for the reader's convenience. Represent in polar 
variables as (j){r,6). Let (f){r,n) denote the Fourier series representation in the 6- 
variable, i.e., 

(24.3) ^{r,n) = [ e-'''^^ir,e)de. 
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The first step is to obtain f l24.4p . the bound for nonzero frequencies. By the Plancherel 
theorem, 

f.27r 

|0(r,^) -0(r, 0)1^6/^ 



2 



= ^|0(r,n)| 

nGZ 

< r \de(P{r,9)\^d9. 
Je=o 

Integrating in rdr from r = to r = D, we obtain the desired bound for nonzero 
angular frequencies: 

(24.4) ||0(r,^) -0(r,O)|U2^(,<z,) < | IV^IU^^ . 

The next step is to obtain (124. 7p . the bound for zero angular frequencies. By the 
mean- value theorem for integrals, there exists tq G [|, 1] such that 

0(ro,O) = 2 [ 0(r,O)rfr. 

Jr=l/2 

By Cauchy-Schwarz 

1/2 



1/2 

By the definition of 0(ro, 0) from (124. 3p with n = and Cauchy-Schwarz, 

1/2 

\kr,0)\<^j \(l>{r,9)\'d9 

Combining we obtain 

(24.5) |0(ro,O)| <( r r \4>{r,9)\\drd9\ ' < ||</.i|Li,(.<i) 

\Jr=l/2Je=Q J 

By the fundamental theorem of calculus, 

i'S=r 

0(r,O)-0(ro,O) = / ds4){s,Q)ds. 

Js=ro 
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Inserting the definition of 0(s, 0), 

/s=r /•2tt 
/ ds(j){s, d)ddds. 
_=ro JB=0 



By Caucliy-Scliwarz, 
(24.6) 

(ps=D /.27r \ 1/2 / ps=D p2-K \ 1/2 

/ / \d,(i){s,e)\hdeds] / / s~Hsde) 
Js=ro Je=o / \Js=ro Je=o J 

<(logD)i/lV0|U2^(.<^). 
Combining fl24.5p and fl24.6p gives 

10(^,0)1 < ||0|U.^(,<i) + (log Z})i/lV0|U|^(.<^) . 

Integrating against rdrdO, 

(24.7) \\krM\\ir<D) < D' (||0|li.^(.<i) + (logD)||V</.||i.^(,<^)^ 

Combining fl24.4p and 024.70 . we obtain the claimed bound. □ 

Now we return to proving (124.20 . It follows easily from (124. ip that u(t) — )■ u* 
weakly in L^(M^). Thus, to establish (124. 2p . it suffices to establish 

(24.8) I Kr = lim J \u{t)\'- 

Let R{t) = e''/^WA(t) for some constant < a < 1, so that R{t) \ as t T just a 
little slower that A(t). We first claim that 

(24.9) / \u{t)\^dxdydz 0. 

We now prove ( 124. 9p . By rescaling and applying Lemma 124.11 

\u{t)\'^dxdydz 

|r-r(t)|2 + |z-2(t)|2<ii(t)2 

|e(f, z)\'^ fj.(f)drdz 

f;2_j_22<e2ii/6(t) 

< / \e{r,z)\'^drdz 

J^2_,_22<e2a/6(t) 

< e'"/'(*) ( [ \e{r, ~z)\^dfdz + \ ! | Ve(f, z)\''dfd?\ 

\/r2+22<l Jf:2_,_52<e2a/6(t) / 

2a/b{t) 

< -^S{t) ^ 0. 
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Let 0(r, z) = 1 when + > 1 but 0(r, z) = for + 2;^ < |. It remains to prove 
that 

lim J (j) ^ ^^^j 7 ^^^j ^ [^(t)^^^^?/^^ = y |'u*|^(ix(i?/(iz. 
To prove this, it suffices to show that 



(24.10) 



' i?(t) y ' ^ ^' y ^ \^ R{t) ' 



^ 0. 



(Note that here we have replaced u by u; this is acceptable, since the contribution 
from Qb is negligible in this region of space.) To prove this, let, for t fixed, 

v{r) = 1 (^^^^^, Hr)\^dxdydz. 

Then observe 

(rr,Zr) f 2 f 

drV = 7^7T~ ' / ■)\'^\^ dxdydz + ——Im / ■) ■ Vu udxdydz. 

R{t) J R{t) J 

By dEZl), \{rr,Zr)\ < X-\ and by BSO 3, \\Vu\\l2 ~ X-\l + 0{T^/^)) ~ A-^. Hence, 

1 1 



\drV{T) 



< 



- e'^imxit) A(r) ■ 



Note that by (I22J]), 

- l\T-T)-y\\og\\og{T-T)\Yl'(l+ \ -) 

it V log|log(T-r)| log(T-r)y 

= 2(T-t)i/2(log|log(T-t)|)V2. 
The rate of h{t) in (|22:8|) implies that e"/^^*) ~ | log(T - t)\^ for some C > 0. Hence, 

' ^^ ^Wl~ga/fe(t) A(t)A A(r) ~ |log(T-t)P ■ 

But by fl24A|) . 

V — r(T) z — z(T) 



f (T) = lim v{t) = d i - 

r^T ' J \ 



R{t) ' R{t) 



yielding ( 124. lOp . and hence, (124. 2p is established. This finishes the proof of Theorem 

o 
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